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ABSTRACT 


\ 

The  purpose  of  this  study  is  to  investigate  and  eval- 
uate adaptive  tracking  filters  for  shipboard  gun  fire  ccmtrol 
systems  that  must  defend  against  evasive  targets.  Several 
design  concepts  are  compared  using  both  tracking  error  and 
predicted  position  error  as  performance  measures.  Tlie 
effects  of  target  evasion,  sensor  measurement  noise  level, 
modeling  uncertainties,  and  length  of  the  measurement  inter- 
val are  investigated,  an ! the  trade-offs  between  performance 
and  algorithm  complexity  are  discussed. 
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1.  INTRODUCTION 

1.1  BACKGROUND:  GUN  HRE  CONTROL  SYSTEMS 

The  need  to  develop  snipboard  gun  fire  control  systems  that  are  capa- 
ble of  destroying  high-speed,  maneuverable  enemy  missiles  creates  a require- 
ment for  pointing  and  tracking  techniques  that  more  accurately  account  for 
target  motion  than  do  conventional  systems.  In  principle,  this  can  be  accom- 
plished by  hardware  improvements,  such  as  more  accurate  tracking  sensors 
and  higher  bandwidth  gun  control  loops,  which  allow  the  gun  pointing  line  to 
respond  more  quickly  to  target  maneuvers.  However,  an  retractive  alternative 
approach  that  does  not  depend  upon  advances  in  hardware  technology  is  the  use 
of  modern  estimation  and  control  theory  to  develop  sensor  data  processing 
techniques  — i.e.,  computer  software  --  which  Uke  maximum  advantage  of 
the  known  mathematical  model  of  target  motion.  With  respect  to  conventional 
we^xm  systems,  this  represents  improved  softu’are,  rather  than  hardware  -- 
realized  by  improving  the  computational  algorithms  used  to  generate  gun 
pointing  commands. 

The  principal  elements  of  a fire  control  system  are  illustrated  in  Fig. 
1.1-1  for  a single  axis.  Briefly,  one  or  more  tracking  sensors  provides  target 
position  measurements  vdiich  are  processed  in  a computer,  together  with 
measurements  of  ship’s  motion,  to  obtain  estimates  of  target  position,  velocity, 
and  acceleration  in  inertial  coordinates.  The  latter  are  appropriately  combined 
to  determine  the  drive  signal  for  bringing  the  actual  pointing  line  mto  coinci- 


dence with  the  commanded  pointing  line. 
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• dynamic  lag  attributed  to  the  finite  bandwidth  of 
the  gun  servo  which  prevents  achieving  exact  equality 
between  the  actual  and  commanded  pointing  lines 
against  a moving  target. 

• The  prediction  error  caused  by  the  fact  that  the 
comm^'^ded  pointing  direction  is  not  exactly  in  the 
direction  required  to  hit  the  target  because  of  target 
trajectory  prediction  errors.  The  latter  are  in  turn 
caused  by  imperiect  modeling  of  target  motion  and 
tracking  sens^ji*  measurement  noise. 

• The  projectile  ballistic  dispersion  error  produced 
by  unknown  aerodynamic  effects  along  the  projectile 
flight  path,  as  w\\  as  non-uniformity  in  projectile 
characteristics  and  firing  conditions. 
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To  some  extent  the  above  errors  will  be  reduced  as  specifications  on 
coniponent  and  instrument  quality  become  progressively  more  stringent.  How- 
ever an  important  alternative  (and  potentially  less  costly)  route  to  system 
imv^»i»vement  is  to  construct  computer  algorithms,  using  estimation  theoiy., 
whi(*h  extract  the  maximum  amoun  of  useful  information  from  the  tracking 
mr  n.surement  data,  thereby  reducing  the  prediction  error.  Tliis  is  especially 
imp.-rtant  in  the  case  of  highly  maneuverable  t’D.rgets  where  prediction  error 
car'  -asily  be  the  most  significant  error  source  in  the  entire  gun  fire  control 
sy;-  The  intent  of  this  report  is  to  investigate  the  potential  for  better  pre- 
di«  on  accuracy  through  software  improvements. 

The  heart  of  the  computer  in  Fig.  1.1-1  is  the  target  tracking  filter 
V • combines  the  tracking  measurements  to  estimate  the  target  position, 
vc*i  H ily,  etc. , with  respect  to  a stabilised  coordinate  frame.  The  latter  are 
1 red  to  ultimately  calculate  the  projectile  time  of  flight,  the  predicted 
tai  ; position  at  the  impact  point,  and  the  gun  pointing  commands.  Conse- 
Qv  .'ly,  it  is  desirable  that  the  tracking  filter  produce  estimates  whose  errors 
: s small  as  possible,  in  some  sense.  Modern  estimation  theory  provides 
a tematic  procedure  lor  accomplishing  this  goal.  Basically,  if  the  target’s 
nr  • • n and  the  tracker  measurement  errors  can  be  described  by  appropriate 
su.  . uciil  mathematical  models,  then  recursive  digital  algorithms  such  as  the 
K “1  filter  are  available  which  will  yield  minimum  variance  estimates  of  the 
y:  ‘ ' les  (called  state  variables)  which  describe  the  target’s  motion.  Such  an 
'4hm  is  called  an  optimal  filter;  its  dependency  upon  the  target’s  equations 
jiiun  and  parameters  and  the  noise  statistics  is  emphasized  in  Fig.  1.1-2. 

The  concept  of  optimal  filtering  is  currently  under  investigation  by  the 
N.r.  . for  both  the  MK  68  and  MK  86  gun  fire  control  systems.  In  future  appli- 
a pix)blem  of  particular  concern  is  the  degradation  in  tracking  filter 
p I jnance  observed  when  a target  performs  significant  unexpected  maneuvers. 
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Figure  1.1-2 


Illustration  of  Factors  Entering  into 
Tracking  Filter  Design 


As  stated  above,  the  filter  construction  requires  knowledge  of  the  mathematical 
model  of  target  motion;  if  the  actual  motion  violates  the  model,  then  the  filter 
is  no  longer  optimal.  Modeling  errors  can  also  arise  with  respect  to  the  tracker 
measurement  noise  statistics.  For  example,  the  errors  in  radar  measurements 
are  to  some  extent  dependent  on  the  target  reflection  properties  which  are  fre- 
quently not  well  known.  Hence,  the  statistical  parameters  (rms  noise  levels) 
needed  for  designing  the  filter  will  be  in  error. 


In  practice,  modeling  errors  are  unavoidable  because  complete  infor- 
mation about  a maneuvering  target  and  the  tracking  sensor  measurement  noise 
is  not  available.  Thus  optimal  filter  design  is  impossible  and  the  problem  is  to 
design  a suboptimal  tracking  filter  which  gives  estimates  that  are  close  to  opti- 
mum. One  common  approach  to  suboptimal  design  is  to  experimentally  select 
the  fixed  filter  which  gives  the  oest  performance  over  the  complete  range  of 
tracking  situations  that  will  be  encountered,  i.e.  it  is  the  best  design  ”on  the 
average*'.  Another  approach  is  to  use  an  adaptive  design  of  the  type  illustrated 
in  Fig.  1.1-3,  in  wiiich  auxiliary  real-time  computations  are  performed  to 


\j 
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Figure  1.1-3  Structure  of  an  Adaptive 
Tracking  Filter 

obtain  better  knowledge  of  the  model  parameters,  and  thereby  improve  the 
filter  estimation  accuracy  on-line.  The  purpose  of  this  study  is  to  investigate 
adaptive  filtering  techniques  that  are  designed  to  track  highly  manuverable 
missile-type  targets,  and  to  compaie  their  performance  with  conventional  fixed 
filters  over  a range  of  target  engagement  situations. 


1.2  TECHNICAL  APPROACH 

In  this  report,  a realistic  but  comparitively  simple  target  tracking 
problem  is  defined.  The  optimum  filter  for  this  problem  is  discussed  and 
several  practical  suboptimal  designs  are  developed  in  detail.  These  include 
both  adaptive  and  fixed  filters  which  were  selected  after  a careful  review  of  the 
technical  literature;  they  represent  original  work  as  well  as  applications  of 
available  techniques.  A comparitive  study  of  these  designs  is  conducted  using 
two  related,  but  distinct,  performance  measuies.  First,  they  are  evaluated  on 
the  basis  of  estimation  error  — i.e.,  the  error  between  the  target’s  slate 
variables  (position,  velocity,  etc.)  and  their  estimates.  This  indicates  how 
well  each  filter  can  "track”  the  current  behavior  of  the  target.  Since  the  ultimate 
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I 

I I task  of  the  tracking  system  is  to  predict  the  future  position  of  the  target  so  that 

i a gun  pointing  line  c m be  calculated,  the  various  designs  are  also  compared 

I I based  upon  their  prediction  error  i.e. , the  error  between  the  future  target 

po  >ition  and  its  predicted  position,  "^he  sensitivity  of  each  design  to  changes 
j in  target  manuver  behavior,  range,  ^ sensor  noise  are  studied  as  well  as  the 

( 

effects  of  noise  correlation  and  measuren*  mt  data  rate.  The  results  indicate 
‘ the  ultimate  prediction  accuracy  possible  with  the  various  tracking  filters  for 

different  target  maneuver  characteristics,  and  the  tradeoff  between  algorithm 
complexity  and  tracking  accuracy. 

) 

\ 

1 . 3 ORGANIZATION  OF  REPORT 

! 

In  Chapter  2,  the  tracking  problem  is  formulated  and  the  assumptions 
and  simplifications  made  for  this  work  are  discussed.  Chapter  3 describes  the 
design  of  tracking  filters  based  upon  optimal  filtering  theory,  discusses  practical 
suboptimal  fixed  configuration  filters,  and  describes  the  use  of  adaptive  tech- 
niques as  a means  of  obtaining  better  tracking  accuracy.  Chapter  4 presents 
the  simulation  results  obtiiined  when  each  design  is  tested  over  a range  of 
tracking  conditions  and  its  performance  is  measured  in  terms  of  estimation 
error  and  prediction  error.  A summary  of  the  results  and  conclusions  are 
presented  in  Chapter  5. 

The  appendices  provide  background  material  for  the  filtering  techni- 
ques investigated.  Appendix  A summarizes  optimal  filtering  and  prediction 
theory  and  Appendix  B discusses  the  details  of  designing  adaptive  filters. 
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TRACKING  PROBLEM  FORMULATION  3 

i 

li 

[ i ! 2.1  TARGET  TRACKING  AND  PREDICTION  j 

r ! 

' i 

. I Li  this  chapt^^r,  the  spo/*tfic  tracking  problem  to  be  studied  Is  devel-  i 

oped  and  discuss. ed.  This  simplified,  but  realistic,  formulcition  is  the  basis 
upon  which  tracking  filters  are  designed  in  Chapter  3 aiid  evaluated  in  Chapter 

4. 

/ 

The  principal  elements  of  a two-dimensional  tracking  problem  are 
I illustrated  in  Fig.  2.1-1.  The  radar  (or  other  tracking  sensor),  for 


Rf16$0 
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; , simplicity  assumed  to  be  fixed  in  inoii:ial  space,  takes  noisy  measurements 

I ^ of  the  target  range  (r)  and  bearing  (B)  at  regular  intervals.  The  tracking  filter 

! ^ processes  these  measurements  and  estimates  the  target’s  state  variables  — 

i f 

j ^ ^ i.e, , those  variables  which  describe  the  target’s  motion  such  as  its  position, 

j velocity,  acceleration,  etc.  Given  these  estimates  the  tracking  filter  is  then 

used  to  predict  the  future  position  of  the  target  by  propagating  the  dynamic 
equations  ci  the  target  forward  in  time  from  the  current  state  estta  js.  The 
total  error  in  prediccod  pt-)sition  depends  upon  the-  length  of  the  prediction  inter- 
val (the  projectile  flight  time),  uiiknown  inputs  io  the  target  dynamics  during 
this  interval,  errors  in  the  target  prediction  model,  and  errors  in  the  current 
^ estim?t  :s  of  the  target’s  states.  The  ultimate  objective  is  to  minimize  the 

prediction  error.  The  length  of  the  prediction  interval  is  determined  by  the 
j target’s  trajectory,  firing  delays,  and  projectile  velocity.  With  the  future 

inputs  to  the  target  dynamics  unknown,  the  only  way  to  reduce  the  prediction 
I error  is  to  model  the  target  dynamics  as  accurately  as  possible  and  design  a 

tracking  filter  which  gives  the  best  possfljle  estimates  of  the  target’s  current 
! state. 

j Modern  optimal  estimation  techniques,  such  as  the  Kalman  filter,  can 

be  applied  to  the  above  tracking  problem.  However,  they  must  be  used  with 
some  care  because  the  filter  design  requires  complete  knowledge  of  the  target 
dynamic  equations  and  the  statistics  of  all  random  inputs.  In  this  problem,  the 
j control  policy  or  attacking  strategy  of  the  target  is  assumed  to  be  partially  or 

completely  unknown.  The  target  might  be  following  a deterministic  guidance 
j law  or  it  might  be  taking  evasive  action  in  a rsindom  manner.  Another  poten- 

tially unknown  quantity  is  the  level  of  measurement  noise,  which  can  depend 
upon  atmospheric  conditions  and  upon  the  size  J»nd  shape  of  the  target,  its  range, 
and  its  changing  reflection  properties  as  it  maneuvers.  With  all  of  these 

I ♦ 

Known  sensor  motion  relative  to  inertial  space  can  easily  be  subtracted  out 
of  the  problem;  hence  no  generality  is  lost  with  this  assumption. 
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uncertain  elements  in  the  problem,  it  is  impossible  to  design  an  optimum 
tracking  filter.  However  modem  estimation  theory  provides  the  basis  upon 
which  suboptimal  filters  can  be  designed  to  give  fairly  good  performance  over 
a wide  range  of  possible  conditions;  the  application  of  this  technology  to  the 
fire  control  tracking  problem  is  the  goal  of  this  report. 

The  remainder  of  thk'  chapter  formulates  in  detail  the  particular 
tracking  problem  to  be  hlmied  i\nd  defines  the  conditions  under  which  the 
tracking  filters  will  be  designed.  The  simplifying  assumptions  used  in  this 
work  are  also  discussed. 


2. 2 CHOICE  OF  TRACKING  COORDINATE  SYSTEM 

In  this  section  the  target^s  equations  of  motion  and  the  tracking 
measurement  equations  are  expressed  in  both  rectangular  and  polar  coor- 
dinates and  a simplified  target  motion  model  based  on  polar  coordinates  is 
discussed. 


Referring  to  Fig.  2.2-1,  the  equations  of  motion  of  the  target  defined 
in  rectangular  coordinates  arc 


X = a 


y = a 


(2.2-1) 


Assumkg  that  the  radar  or  other  tracking  svmsor  takes  noisy  measuremenis, 
and  z^,  of  the  target's  range  and  bearing,  the  measurement  equations  are 


2 2 

Zj.  = fx  +y  ^v^ 
Zg  = tan"*  (y/x)  + 


(2. 2-2) 
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Figure  2.2-1  Target  Acceleration  Component  Definitions 


where  v and  v are  measurement  errors.  The  range  measurement,  z , is 
r 6 r 

the  true  range,  r,  plus  the  range  measurement  noise,  v , and  z is  similarly 

r u 

related  to  the  true  bearing,  0.  Note  that  while  the  dynamic  equations. 

Eg.  (2.2-1),  are  linear  in  the  rectangular  (x,  y)  coordinates,  the  correspond- 
ing measurement  equations,  Eg,  (2.2-2),  are  highly  nonlinear. 

To  derive  the  equivalent  set  of  equations  in  polar  coordinates  the 
following  transformation  is  used; 


X = r cos  0; 


y = r sin  0 


(2.2-3) 
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Differentiating  each  of  these  equations  twice  with  respect  to  time  and  solFing 
the  resulting  pair  of  simultaneous  equations  for  r and  e gives 

r = re^  + ic  cos  0 + y sin  3 

(2.2-4) 

••  -2fe  . -X  sin  6 + y cos  6 
r r 

Note  that  the  target's  acceleration  components  along  and  perpendicular  to  the 
Une-of-sigjit  (see  Ffg.  2.2-1)  are  respectively  given  by 

a^  = X cos  0 -I-  y sin  6 

{2.2-5> 

- -X  sin  e + y cos  & 


Substitution  from  Eq.  (2. 2-5)  into  Eq.  (2. 2-4)  produces 

••  *2 

r = re  + 


(2.2-6) 


The  measurement  equations  are 


z = r + V 
r r 


(2.2-7) 


(2.2-8) 


Thus,  in  polar  (r,9)  coordinates  the  syfitem  has  nonlinear  dynamic  equations 
linear  measurement  eqi  it  ions . 


For  the  purpose  of  evaluating  different  target  tracking  filters,  a sim- 
plified tracking  model,  defined  in  terms  of  polar  coordinates,  is  employed. 
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With  respect  to  the  target's  angular  motion,  ithe  tenn  involving  r in 
Eq.  (2. 2-7)  is  neglected,  yielding 


l[2.2-9) 


This  assumption  is  justified  for  evalmting  design  concepts  because  the  pri- 
mary error  in  tracking  6 tends  tc  be  that  caused  by  a , the  unknown  tairget 

B 

normal  acceleration  — i.e.,  its  acceleration  perpendicubtr  to  the  line-of- 
sig^t.  Likewise  in  Eq.  (2. 2-6)  the  primary  error  in  tracking  range  Is  caused 
by  a^,  and  this  equation  can  be  approximated  by 


r = a 


(2. 2-10) 


Since  Eqs,  (2. 2-9)  and  (2. 2-10)  have  the  same  form  it  is  sufficient  to  test 
adaptive  tracking  schemes  using  Eq.  (2. 2-9)  as  the  simplified  model  off  the 
true  system.  Once  a promising  adaptive  scheme  has  been  found  then  it  can  be 
easily  applied  to  tracking  a target  that  actually  obeys  Eqs.  (2. 2-6)  and  (2. 2-7). 

In  summary,  the  various  filters  in  this  report  are  tested  against  a 
target  that  obeys  Eq.  (2. 2-9).  These  tests  are  made  at  a constant  range,  r, 
and  sensitivity  studies  are  performed  to  determine  how  the  estimation  errors 
vary  with  range. 


2.3  MEASUREMENT  NOISE  AND  TARGET  REFERENCE  MODEL 

To  test  tracking  filters  using  a computer  simulation,  a measurement 
noise  and  target  reference  model  must  be  defined.  In  an  actiial  tracking 
situation  the  measurement  data  comes  from  the  sensor  and  only  the  filter  is 
implemented  in  the  tracking  computer.  In  a computer  simulation,  the  source 
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of  the  measurement  data  must  be  simulated  as  well.  The  measurement  noise 
and  target  reference  model  performs  this  function.  In  this  report  the  target 
reference  model  consists  of  the  simplified  tracking  model  Eq.  (2.2-9),  and  a 
variety  of  profiles  for  the  target’s  normal  acceleration,  a . The  measurement 

u 

noise  reference  model  consists  of  errors  calculated  using  the  computer  random 
number  generator.  These  errors  are  then  added  to  the  target's  computed 
angular  position,  e,  to  obtain  noisy  measurements,  z . 

Figure  2.3-1  is  a block  diagram  of  the  target  and  measurement 
reference  model.  For  a particular  test  run  the  following  must  be  defined: 


r - target  range,  (ft) 

• 3 

a (t)  - normal  acceleration  rate  time  history,  (ft /sec  ) 

s 

At  - measurement  interval,  (see) 

V (t  ) - measurement  noise  sequence,  (rad) 

s n 

6(0)  - initial  angular  position,  (rad) 

t^(0)  - initial  angular  velocity,  (rad/sec) 

2 

a (0)  - initial  normal  acceleration,  (ic/sec  ) 

6 


Figure  2.3-1  Target  and  Measurement  Reference  Model 
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The  type  of  tarp;et  trajectory  produced  by  this  model  is  determined  by 

the  choice  of  a^Ct).  Figure  2.3-2  gives  three  examples  of  acceleration  rate 

profiles  that  might  be  produced  by  a target  with  a large  maneuver  capability. 

Both  a and  the  resulting  a are  plotted  for  a ten-second  trajectory.  Part  (a) 

6 6 

of  this  figure,  where  the  target  has  a constant  normal  acceleration,  represents 

what  will  be  called  the  nonevasivr  case.  Part  (b)  is  referred  to  as  a ’’mildly” 

evasive  target  which  changes  its  acceleration  with  the  maximum  absolute  value 

3 3 

of  a being  50  fl/sec  and  the  rms  value  being  30  ft/sec  or  about  1 g/sec.  This 

is  called  the  nominal  case.  Finally,  i^irt  (c)  of  the  figure  is  the  highly  evasive 

case  which  represents  tho  worst  evasive  motion  considered  in  f study  and 

anticipates  the  possible  capabilities  of  future  antishippLng  missile  threats. 

3 

Here  the  maximum  absolute  \^lue  of  a is  150  ft /sec  and  the  rms  value  is 
3 ® 

90  ft/sec  or  approximately  3 g/sec.  The  filters  in  this  report  arc  optimized 
for  the  nominal  case,  and  then  tested  for  all  three  cases  to  determine  their 
performance  over  the  complete  range  of  target  behavior. 

Choosing  the  noise  sequence,  v (t  ),  defines  the  reference  model  for 

S n 

the  measurement  noise  process.  The  set  of  times,  l^,  are  the  sampling  or 
measurement  times.  Two  types  of  models  are  used  in  this  report.  The  first 
is  an  uncorrelated,  or  \ndiite,  gaussian  sequence.  Here  the  samples  are 
cliosen  from  a zero  mean  gaussian  random  number  generator  where  the  indi- 
vidual members  of  the  sequence  are  indej>endent  of  each  other,  and  the  rms 
value  of  the  sequence  is  That  is. 


E|v,(t„)) 


for  all  n 


n / m 


hy  n = m 
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rmt  vjlu« 


TlMt  (iecJ 


(ai  Noari/asive  Case,  Constant  Acceleration  (a 


tmt  value  • 30  fi/iec' 


TIME  (mc) 


150  ft/scc3) 


Ic)  Higtily  FvJSfvt?  C*isc  ( !a, 


Fir^urc  2.3-2 


Tarj^et  Norma'  Acceleration  Profiles  and 
Their  Derivatives 
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The  second  noise  model  is  a correlated  or  ’’colored”  paussian  sequence,  having 
correlated  samples  generated  by  passing  a while  sequence  like  that  defined 
above  through  a lov/  pass  filter.  Specifically,  the  measurement  noise  sequence 
is  generated  by  driving  the  difference  equation  for  a low-pass  filter  with  the 
output  of  a gaussian  random  number  generator.  The  equation  is 


(2.3-1) 


where  is  a zero  mean  gaussian  white  noise  sequence,  r is  the  correlation 

time,  At  is  the  measurement  interval,  and  v (t  ) is  the  measurement  noise  at 

d n 

time  t^.  To  keep  the  rms  value  of  v^  (t^)  equal  to  the  nominal  value,  the 
rms  value  of  w^,  denoted  c,  is  given  by 


a = t/ 1 - e 


-2  At  A 


(2.3-2) 


Thus  V (t  ) is  a correlated,  or  ’’colored”,  gaussian  noise  sequence  with  an 
tf  n 

rms  value  equal  to  and  a correlation  time  of  seconds. 

In  the  tracking  problem,  the  sensor  (e.g. , radar)  receiver  noise  is 
usually  wide-band  so  that  white  gaussiaii  noise  is  a good  model.  However  the 
total  measurement  error  can  also  contain  low  frequency,  correlated  noise 
components  caused  by  a slowly  varying  target  cross  section,  due  to  the  target’s 
own  motion.  These  considerations  motivate  the  two  measurement  noise  models 
above. 
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3.  TIUCKIWG  FILTER  DESIGN 

3.1  OPTIMAL  FILTER  DESIGN  AND  PERFORMANCE 

The  best  tracking  algorithm  design  possible  in  terms  of  achieving  the 
minimum  mean  square  estimation  error  is  the  Kalman  filter  outlined  in 
Appendix  A.  This  design  requires  that  the  tracking  system  dynamic  and 
measurement  equatic  s be  completely  known  and  linear.  In  addition,  the 
driving  functions  to  the  system  differential  equations  must  be  white  gaussian 
noise  processes  with  known  mean  and  spectral  density  and  the  measurement 
noise  must  be  a gaussian  sequence  with  kno\vn  mean  and  jo variance.  If  these 
conditions  are  met,  the  tracking  filler  design  is  completely  determined. 

In  this  study,  the  filter  designs  are  based  upon  the  simplified  linear 

model 


e = -jr  (3.1-1) 

^ (3.1-2) 

Where  the  range  r is  assumed  known.  Equations  (3.1-1)  and  (3.1-2)  have 

the  same  form  as  the  reference  model  defined  in  Chapter  2;  however  the  input 

to  the  reference  model  — the  target  acceleration  a in  Fig.  2.3-2  --  is  not 

accurately  represented  as  while  noise  and  the  statistical  properties  of  the 

measurement  noise  sequence,  v (t  ),  may  not  be  kno\vn  accurately.  These 

6 n 

facts  prevent  (lie  desipi  of  an  optimal  tracking  filler. 
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It  is  possible  to  artificially  design  an  opthnum  filter  the  reference 
model  by  assuming  that  the  target  acceleration  rate  is  known  and  the  measure- 
ment noise  is  gaussian  white  noise  with  known  covariance.  An  analysis  of  this 
case  gives  «.  measure  of  filter  performance  under  ideal  conditions.  Equations 
(3.1-1)  and  (3.1-2)  are  the  system  dynamic  and  measurement  equations,  where 

a and  r are  given  and  the  covariance  (in  this  case  the  mean  square  value)  of 
® 2 

v^  (tn)  is  The  dynamic  equation  can  be  put  into  state  variable  form  cor- 
responding to  Eq.  (A.  1-1)  of  Appendix  A by  defining 


’Xj(t)- 

-o(t)  ■ 

x(t)  = 

X2(t) 

= 

e(t) 

-01  0 ■ 

"o’ 

F(t)  = F = 

0 0 1/r 

G(t)u(t)  = 

0 

0 C 0 

.0. 

and 


0 


0 

a^(t) 


Then  Eq.  (3.1-1)  becomes 


x(t)  = Fx(t)  f b (t) 


(3.1-3) 
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The  measurement  equation  is  put  in  the  form  of  Eq.  (A.  1-8)  by  defining 


Then 


z = 
n 


z^(t  ) 
0 ' n 


Vn  = 


H = = [1  0 0] 


Z 

n 


Hx(tn)  + V 


n 


(3.1-4) 


Since  measurements  are  taken  only  at  discrete  instants  of  time  it  is  convenient 
to  express  Eq.  (3.1-3)  in  discrete  time  also.  It  is  assumed  that  the  tjacking 
filter  is  implemented  on  a digital  computer  so  this  formulation  will  be  most 
natural  and  efficient. 


Assuming  a uniform  measurement  interval  of  length  At  and  using  the 
definitions  of  Appendix  A,  Eqs.  (3.1-3)  and  (3.1-4)  become 


X = <^x  - + b ^ 

-n  ^-n-1  -n-i 


z 

n 


+ V 


n 


(3.1-5) 


where 


4>  = 


FAt 

e 


t -At 
n 


F(t„-r) 


b (r)dT 
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To  complete  the  definitions  required  for  the  filter  equations  note  that  the  driv- 
ing noise  covariance  matrix.  0^  is  equal  to  zero  (see  Eq.  A.  1-7)  and 


R = E[v  ] = 
n ^ n^ 


(1x1  matrix) 


where  c/^  is  the  rms  measurement  noise  level  of  the  reference  model.  By 
assuming  values  for  the  filter  initial  conditions,  Xq(-)  and  Pq(-),  the  filter 
equations  given  in  Fig.  A.  1-1  are  completely  defined.  Again,  it  is  empha- 
sized that  both  a as  a function  of  time  and  a , the  rms  value  of  measurement 
0 V 

noise  used  in  the  reference  model,  must  be  known  to  implement  the  optimum 
filter. 


The  filteiing  algorithm  for  the  case  defined  above  is 


X 

-n 


= ^ ^n-1  ^ K ' n-1  + ^n-l^  ^ 


(3.1-6) 


where  is  the  optimum  state  estimate  of  immediately  after  the  measure- 
ment is  taken  and  k^  is  the  Kalman  gain  vector  defined  by  Eqs,  (A.  1-14), 
(A.  1-15)  and  (A.  1-16).  The  filter  performance  is  shown  in  the  figures  that 
follow  and  is  based  upon  the  following  parameter  values: 


1 . 4 mrad 
6.1  sec 
10, 000  ft 
0 


(rms  measurement  noise  level) 
(measurement  sampling  interval) 
(target  range) 

(initial  state  estimate) 


(3.1-7) 


! 

I 
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with 


'2.0x  10"^ 

rad 

x(0)  = 

4.0  X 10"^ 

rad/sec 

. 2.0  X 10^ 

ft/sec^  . 

2.0  X 10”®rad^ 

0 

0 7.1 

X 10"^  (rad/sec)^ 

0 

0 

2.7  X 

(3.1-7) 

cont. 

0 

0 

10^(ft/sec^)^ 


Figure  3.1-1  shows  the  target  position,  velocity,  and  acceleration  esti- 
mation errors  as  functions  of  time  for  a ten-second  simulation  using  the 
nominal  case  for  a^  shown  in  Fig.  2.3-2(b).  Part  (a)  of  the  figure  shows  the 
angular  estimation  error  on  the  right  vertical  scale  and  the  position  estimation 
error  (measured  normal  to  the  line-of-sight)  on  the  left  vertical  scale.  For  a 

range  of  10, 000  ft.  the  small  angle  approximation  e.  = sin  (e^)  is  valid  since 
^ 0 6 
e = e - 6 is  much  less  than  5 degrees.  Thus 


1 

3i 

I 


re^  = r sin  (e^) 

gives  the  position  error  in  feet.  Part  (b)  of  the  figure  makes  the  same  approx- 
imation lo  obtain  the  velocity  error  from  the  angular  velocity  error. 

Since  the  target  acceleration  rate  is  assumed  known,  it  is  accounted 

for  exactly  by  the  term  b^_^  in  the  filter  equation,  Eq.  (3,1-6);  therefore  the 

errors  shown  in  the  figure  are  independent  of  the  particular  choice  of  a (t). 

0 

This  is  a basic  property  of  the  optlinai  filter  since  the  equtions  which  propa- 
gate and  update  the  estimate  of  the  error  covariance  matrix,  P^,  (Eqs.  (A.  1-15) 
and  (A.  1-16)),  are  independent  of  the  system  input  vector  b(t),  which  contains 
a^(t).  Part  (d)  of  Fig.  3.1-1  shows  how  the  mean  square  estimation  errors 
converge  towards  /ere  as  calculated  by  the  filter  covariance  equations.  These 
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curves  are  in  qualitative  agreement  with  the  actual  estimation  errors  shown  in 
parts  (a),  (b),  and  (c)  of  the  same  figure. 

Note  that  the  estimation  errors  in  Fig.  3.1-1  parts  (b)  and  (c)  become 
small  relative  to  their  initial  values  after  approximately  two  seconds.  The 
rms  values  given  for  the  errors  are  calculated  using  only  the  last  eight  seconds 
of  the  simulation  to  remove  the  effect  of  this  initial  transient.  For  the  purposes  of 
target  prediction  discussed  later,  notice  that  the  three  estimation  errors  are 
highly  correlated;  that  is,  when  the  normal  acceleration  estimation  error  is 
large  and  positive,  so  are  the  position  and  velocity  estimation  errors.  This 
tends  to  cause  the  components  of  the  rms  prediction  error,  contributed  by  the 
individual  estimation  errors,  to  combine  additively.  Section  4.3  develops  this 
point  in  more  detail. 

Figure  3.1-2  gives  a clearer  illustration  of  the  convergence  of  the 
normal  acceleration  estimate  as  a fimction  of  time  relative  to  the  actual  value 
of  normal  acceleration.  Note  the  rapid  convergence  to  a small  percentage 
error. 


Figure  3.1-2  Normal  Acceleration  and  its  Estimate  for  Optimum 
Filter:  Nominal  Trajectory 
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The  resultj  in  this  section  demonstrate  the  accuracy  with  which  the 
optimum  filter  can  track  the  target’s  states  if  sufficient  knowledge  is  available 
to  construct  the  filter.  In  the  next  section,  practical  designs  will  be  discussed 
where  this  knowledge  is  not  available. 


3. 2 PRAC'nCAL  FILTER  DESIGN  AND  PERFORMANCE 

The  last  section  demonstrated  how  an  optimum  filter  would  work  if 

the  target  and  measurement  model  were  accurately  known.  The  primary 

unknowns  which  make  the  optimum  design  impractical  are  a (t),  the  time 

6 

history  of  acceleration  rate,  and  the  true  rms  measurement  noise  level. 
This  section  assumes  that  is  known  and  the  major  problem  is  dealing  with 
the  unlcnown  target  normal  acceleration.  The  problem  of  being  different 
than  its  assumed  vah  e is  left  until  Section  4. 2 where  measurement  noise 
sensitivity  is  discussed. 


One  common  method  of  dealmg  with  the  unknown  input  a is  to  assume 

u 

that  it  can  be  approximated  by  white  noise  and  construct  a Kalman-type  track- 
ing; filter  as  described  in  Appendix  A.  Recall  that  the  reference  model  equa- 
tion is 

X = Fx  + b(t) 


(3.2-1) 


A propo;j3ed  filter  design  is  based  upon  a model  of  the  form 


x_  - Fx  ‘fu(t) 
-m  -m  ' 


(3.2-2) 


where  the  notation  denotes  the  model  representation  for  x(t)  and 


u(t)  = 


0 • 
0 
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The  quantity  u^(t)  is  assumed  to  be  a gaussian  wnico  noise  process  with 

E [Uga)  Ug(T)]  = q6(t-T)  (3  2-3) 

That  is,  the  unknown  input  a^is  modeled  as  gaussian  white  noise  with  a spec- 
tral density  equal  to  q,  (see  Appendix  A).  Recalling  Fig.  2.3-2,  it  is  clear 
that  the  a^  signal  shown  does  not  look  like  wiiite  noise.  However,  the  /hite 
noise  moJel  has  two  advantages.  First  it  does  not  increase  the  order  (number 
of  states)  of  the  filter  beyond  that  required  for  the  optimum  filter  discussed  in 
the  last  section.  Second,  if  the  target  has  a constant  acceleration,  then  b(t)  = 
0 in  Eq.  (3. 2-1)  and  can  be  modeled  exactly  in  Eqs.  (3. 2-2)  and  (3. 2-3)  by 
setting  q = 0.  Furthermore,  the  white  noise  model  is  a good  approximation  to 
a^  ii  the  effective  bandwddth  of  the  latter  is  large  relative  to  other  dynamic 
effects. 


The  filter  design  based  on  Eqs.  (3.2-2)  and  (3.2-3)  is  the  same  as 
that  discussed  in  Section  3. 1 except  that  b(t)  is  unknown  and  not  included  in  the 
filter  equations  and  the  matrix  Q^,  which  originally  was  zero,  now  becomes 


0 0 
0 0 
0 q 


\e  / dr 


The  filter  algorithm  is 


+ k fz  - H*x  ,] 

— n '■  n — n-l^ 


(3.2-4) 


where  k^  is  the  Kalman  gain  sequence  defined  by  Eqs.  (A.  1-14),  (A.  1-15),  and 
(A.  1-16)  when  they  are  based  upon  the  design  model.  Note  that  the  Kalman 
gain  sequence  depends  upon  the  choice  of  q.  This  algorithm  will  be  called  the 
suboptimal  third-order  fixed  filter. 

The  value  for  the  spectral  density,  q,  chosen  for  the  suboptimum  filter 
in  Eq.  (3. 2-4)  can  be  any  positive  number;  however  the  particular  value  chosen 
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effectively  determines  the  filter  bandwidth  and  is  therefore  very  important  * 

to  the  performance  of  the  filter.  A low  value  of  q results  in  a low-band- 

width  tracking  filter  which  can  accurately  track  targets  that  manuver  very 

little  — i.e. , the  rms  value  of  a is  small  and  target  motion  is  low  frequency 

Q 

in  nature  so  a low-bandwidth  filter  can  be  used  to  eliminate  the  measurement 
noise  from  the  data.  On  the  othe’’  hand,  if  the  target  is  nighly  evasive  and  s 
motion  has  high  frequency  components,  then  a large  value  of  q is  needed  so  that 
the  filter  has  a high  bandwidth  and  can  respond  rapidly  to  the  target  motion. 

In  this  case,  for  a fixed  measurement  rate  and  measurement  noise  level,  even 
the  best  choice  of  q will  give  poorer  tracking  accuracy  than  the  case  where  q 
is  small  and  the  target  manuvers  very  little.  Thus  the  proper  choice  of  q,  and 
consequently  the  filter  bandwidth,  is  dependent  upon  how  evasively  the  target 
maneuvers.  Hereafter,  q will  bo  referred  to  as  the  design  parameter  for  the 
third-order  fixed  filter . 

To  be  more  specific,  Fig.  3. 2-1  shows  the  results  of  a series  of 

simulations  for  the  same  conditions  given  in  Eq.  (3.1-7),  except  that  the  sub- 

optimal  filter  derived  above  is  employed.  The  trajectory  for  a (t)  chosen  here 

6 

is  that  of  Fig.  2.3-2(b),  the  nominal  case.  The  figure  shows  the  mean  square 

estimation  errors,  again  calculated  by  time -averaging  over  the  last  eight 

seconds  of  each  run  to  eliminate  the  effect  of  the  initial  transient.  The  mean 

square  estimation  errors  of  angular  position,  angular  velocity,  and  normal 

.acceleration  are  plotted  for  a wide  range  of  choices  of  q,  the  model  spectral 

density  which  must  be  selected  by  the  filter  designer.  All  three  curves  are 

minimized  at  approximately  the  same  value  of  q.  Thus,  there  is  a ’’best” 

choice  for  q for  this  trajectory;  but  it  caimot  be  determined  unless  a (t)  is 

0 

known  a priori. 

Figure  3.2-2  shows  how  the  me;in  square  normal  acceleration  esti- 
mation error  varies  as  a function  Ok  q,  for  the  different  target  trajectories 
showai  in  Fig.  2.3-2.  The  bottom  curve  in  Fig.  3.2-2  is  the  nonevasive  case 
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1 10  10'  10^  10^  10® 
DESIGN  PARAMETER,  q |ft^/$cc®) 


Figure  3. 2^-2  Trajectory  Sensitivity  of  the  Third-Order 
Fixed  Filter 

(Fig.  2.3  -2(a)  where  the  target  acceleration  is  constant.  This  curve  indicates, 

as  stated  earlier,  that  the  optimum  choice  of  q is  zero  for  the  nonevasive  case. 

The  middle  curve  corresponds  to  the  nominal  traiectory  (Fig.  2.3-2(b))  and  the- 
3 

best  q is  near  10  . Finally,  the  top  curve  corresponds  to  the  worst -case 
evasive  irajcclory  (Fif;.  2. 3-2(c))  with  the  best  value  of  q ~ 10^,  Recall  that 
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the  optimum  filter  of  Section  3.1  gives  the  same  rms  estimation  errors  re- 
gardless of  the  time-history  of  a^(t)  because  a (t)  is  assumed  known;  this  is 
evidently  not  the  case  for  the  suboptimal  filter  evaluated  in  Fig.  3.2-2.  The 

filter's  performance  deviates  farther  from  optimum  as  the  rms  value  of  a (t) 

0 

increases,  even  if  the  best  choice  of  q is  knowm. 

To  further  clarify  the  operation  of  the  third-order  fixed  filter  design, 
it  is  useful  to  plot  the  target  normal  acceleration  and  its  estimate  as  a function 
of  time  for  the  nominal  trajectory  for  three  different  values  of  q.  Tins  is  done 
in  Figure  3.2-3  for  q equal  io  7.5,  750,  and  75,000  ft  ^sec  . The  rms  esti- 
mation error  differs  significantly  between  these  three  plots.  Part  (a),  which 
is  the  low  q case,  corresponds  to  a filter  having  a low  bandwidth.  Observe 

A 

that  the  estimate  a^  has  difficulty  tracking  rapid  changes  in  a^,  however,  the 
time  history  of  a^  is  very  smooth  indicating  that  most  of  the  measux'ement 
noise  is  suppressed  by  the  filter.  By  comparison.  Fig.  3.2-3(c),  corresponds 
to  a relatively  high  bandwidth  filter.  Here  a^  can  track  a^  when  it  varies 
rapidly,  but  a lot  of  the  measurement  noise  gets  through  the  filter  and  degrades 
the  estimate.  Fig.  3. 2-3(b)  is  a medium  bandwidth  filter  and  represents  a 
nearly  optimum  tradeoff  between  tracking  ability  and  the  amount  of  measurement 

A 

noise  that  corrupts  a^:  thus  it  has  a lower  rms  estimation  error  than  the  other 
two  cases. 


The  third-order  fixed  filter  models  a^  as  white  noise  and  consequently 
as  the  output  of  an  integrator  driven  by  white  noise;  in  this  model  the  target 
acceleration  is  called  a random  w^lk  process.  In  an  effort  to  improve  the 
filter’s  performance,  other  models  for  a^  can  be  investigated.  For  example, 
in  another  set  of  simulations,  a^  \viis  modeled  as  the  output  of  a low  pass  filter 
with  time  constant,  t.  driven  by  white  noise  with  spectral  density,  q;  i.e., 


4 

0 T 0 


(3.2-5) 
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RVSinMOR  b'jOli 


7 4 6 8 

TIME  (trc) 

(a)  Low  Bandwidth  Filtef,  q = 7.5  ft2/sec^ 


300 


(b)  Medium  Bandwidtit  Filter,  q = 750  ftwsec^ 
(Nearly  optin.urr  choice  of  q) 


300 


L-  I »•  M Tl* 


(c)  BuiKlvvidth  ^dtrr.  q 75.000  ll^/vcc^ 


Fi^rc  3.  2-3 


Third-Ordei  Fiiter  NoriDiil  Acceleration  Ivjtimateb 
Nominal  Trajectory 
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In  this  case  two  parameters  must  be  selected,  r and  q.  A filter  based  upon 
this  model  was  simulated  using  the  highly  evasive  trajectory  in  Fig.  2.3-2(c). 
The  mean  square  normal  acceleration  estimation  errors  as  fimctions  of  the 
parameters  in  Eq.  (3.2-5)  are  shown  in  Fig.  3.2-4.  For  comparison,  the 
results  for  the  random  walk  acceleration  model  are  also  given.  It  is  clear 
from  this  figure  that  the  random  walk  model  yields  better  performance  than 
the  more  complex  low-pass  model  just  described.  Note  also  that  the  various 
low-pass  models  are  slightly  more  sensitive  to  the  incorrect  choice  of  q and 
approach  the  random  walk  model  as  r is  increased. 

Another  possibility  is  to  try  a higher-order  acceleration  model.  For 
example,  in  Fig.  3.2-5,  a is  represented  as  the  output  of  a double  integrator 

CF 

driven  by  gaussian  white  noise  with  spectral  density,  q.  Using  this  model,  the 


Figure  3.2-4  Performance  of  Low  Pass  Acceleration  Models: 
Highly  Evasive  Trajectory 
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GAUSSiAN  WHITE 
NOISE  PROCESS 
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SPECTRAL  DENSITY  = ^ 
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Figure  3, 2-5  Doubly  Integrated  Gaussian  White  Noise  Model  for 
Target  Acceleration 


equations  of  Appendix  A give  a fourth-order  fixed  filter  design.  Figure  3. 2-6 
shows  the  rms  normal  acceleration  estimation  errors  achieved  with  this  model 
for  the  three  test  trajectories  and  compares  these  errors  with  those  obtained 
using  the  third-order  fixed  filter  that  assumed  a random  walk  acceleration 
model.  Again  it  is  clear  that  in  each  case,  the  simpler  third-order  filter 
works  better  than  the  fourth-order  filter. 


r 


I. 

i. 


DESIGN  PARAMETER 


{q  (ll»/«c*) 


ORDER  FILTER 
4’‘'OROE  ! JLTER 


Figure  3.2-0  Comparitive  Performance  of  Third-Order  and  Fourth- 
Order  Fixed  Filters  for  Three  Trajectories 
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These  studies  indicate  that  the  fixed  filter  design,  based  upon  the 
random  walk  model  for  a^  shown  in  Fig.  2.3-1  is  preferable  if  the  tracking 
situation  is  similar  to  that  tested  here;  i.e. , for  the  given  data  rate,  measure- 
ment noise  levels  and  target  maneuver  characteristics,  'fhe  main  value  of 
these  studies  is  that  they  show  the  type  of  procedure  that  must  be  followed  to 
design  a filter  for  a particular  class  of  tracking  situations.  Clearly  the  third- 
order  suboptimal  filter  cannot  perform  nearly  as  well  as  the  ideal,  but  im- 
practical optimum  filter  described  in  Section  3.1;  however  it  appears  to  do 
about  as  well  as  possible  if  the  target  maneuver  time-history  is  not  known  a 
priori.  The  primary  design  problem  is  the  choice  of  q since  its  best  value  is 
related  to  the  unknown  rms  value  of  a^  over  the  trajectory.  In  the  next  section 
adaptive  designs  are  considered  in  an  attempt  to  overcome  this  problem. 


3.3  ADAPTIVE  TRACKING  FILTER  DESIGN 

Two  steps  are  important  in  the  design  of  a target  tracking  algorithm. 
First,  a mathematical  model  for  target  acceleration  must  be  chosen  based 
upon  theoretical  considerations  and  experimental  work,  as  described  in  pre- 
vious sections.  This  model  must  be  general  enough  to  cover  all  anticipated 
target  acceleration  profiles.  Second,  a tracking  filter  must  be  designed  based 
upon  this  model,  the  target’s  equations  of  motion,  and  the  measurement  error 
statistics.  Its  performance  should  be  as  good  as  possible  despite  changes  in 
the  target’s  acceleration  profile,  errors  in  the  particular  dynamic  model 
selected  for  target  acceleration,  and  variations  in  the  measurement  noise  level. 
One  approach  is  to  experimentally  determine  the  values  of  the  filter  design 
parameters  that  give  the  best  average  performance  over  all  anticipated  track- 
ing situations:  alternatively  the  parameters  could  be  optimized  for  the  worst 
case — i.e.,  for  the  most  violent  target  maneuvers  and  highest  measurement 
noise  level.  However,  a fixed  parameter  design  based  on  either  of  these 
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procedures  is  likely  to  perform  significantly  worse  than  a design  matched  to 
the  actual  target  maneuver  characteristics  encountered. 

A potentially  better  design  than  the  fixed  parameter  filter  is  an  adap- 
tive filter  which  can  identify  each  particular  tracking  situation  as  it  arises 
and  adjust  itself  for  achievmg  the  best  performance.  For  example,  if  the  type 
of  target  trajectory  could  be  identified  on-line,  then  the  proper  choice  of  the 
bandwidth  parameter  q for  the  third-order  filter  could  be  made.  Many  different 
adaptive  filtering  techniques  have  been  investigated  in  the  past  ten  years  and  an 
extensive  literature  review  was  conducted  to  determine  which  methods  might  be 
useful  for  the  target  tracking  problem  treated  here.  Three  techniques  which 
seem  promising  are  described  in  this  section  and  computer  simulation  results 
are  presented  in  the  next  chapter.  All  of  these  methods  start  by  assuming 
that  the  filter  is  to  have  the  same  basic  structure  as  the  optimal  Kalman  filter 
outlined  in  Appendix  A.  They  differ  in  how  they  change  or  augment  this  1 

structure  to  make  the  filter  adjust  for  modeling  errors  and  changing  conditions. 
References  1 through  17  contain  the  adaptive  filtering  theory  found  to  be  per- 
tinent to  the  tracking  problem.  The  techniques  described  in  this  report  were 
developed  from  this  background  and  represent  both  applications  and  extensions 
of  existing  techniques. 


3.3.1  Adaptive  Bandwidth  Filter 


The  general  structure  of  the  adaptive  bandwidth  filter  is  shown  in 
Fig.  3.3-1,  consisting  of  i\vo  connected  filters.  Filter  1 is  a fixed  configura- 
tion fourth-order  Kajman  filter  which  produces  estimates  of  the  target  evasive 
maneuver  level,  that  is,  the  rate  of  change  of  target  angular  acceleration,  a^. 
Filter  2 is  an  adaptive  band\vidth,  third-order  Kalman  filter  which  estimates 
the  target’s  angular  position,  angular  velocity,  and  normal  acceleration. 

The  bandwidth  of  the  latter  is  controlled  by  the  estimate  of  target  evasive 
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maneuver  level  obtained  from  Filter  1.  Thus  if  the  target  is  not  evading  and 
1 ; is  small,  then  the  estimate  a^  will  be  small  and  Filter  2 is  adjusted  to 

have  a low  bandwidth  to  achieve  good  noise  suppression.  If  the  target  begins 
1 . to  change  its  angular  acceleration  rapidly,  Filter  1 will  reflect  this  fact  in  its 

j estimate  of  acceleration  rate  and  will  accordingly  raise  the  bandwidth  of  Filter 

» . 2.  Essentially  this  adaptive  filter  tries  to  detect  how  much  the  target  is  man- 

^ euvering  and  then  adjusts  itself  to  the  proper  bandwidth  for  that  maneuver 

I level. 

i The  design  of  the  fourth-order  fixed  filter  is  the  same  as  that  discus 

sed  in  Section  3. 2 where  target  acceleration  is  modeled  as  the  output  of  a 
I double  integrator  driven  by  white  noise  with  a spectral  density  designated  by 

the  parameter  q;  see  Fig.  3. 2-5  for  a block  diagram  of  this  model. 

! . 

Recalling  that  the  third-order  filter,  investigated  in  Section  3.2, 

) develops  appreciable  acceleration  estimation  errors  when  the  target  maneu- 

vers evasively,  it  is  expected  that  the  estimates  of  acceleration  rate  from 

). 


FILTER  2 


I 

Figure  3.3-1  Adaptive  Bandwidth  Filter 
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the  fourth-order  filter  will  not  be  very  accurate.  However,  the  philosophy  of 
this  design  does  not  require  accurate  estimates  of  a^.  What  is  needed  is  a 
general  idea  of  its  magnitude  so  that  the  bandwidth  of  the  third-order  filter 
can  be  appropriately  adjusted.  Assuming  that  a rough  estimate  is  available 
for  a^,  the  desiem  of  the  variable  bandwidth  ifilter  in  Fig.  3.3-1  is  now 
discussed. 

Filter  2 is  a third-order  Kalman  filter  which  models  target  accelera- 
tion as  the  output  of  an  integrator  driven  by  gaussian  white  noise  with  a spec- 
tral density  of  q.  For  a constant  q,  this  is  the  same  third -order  filter  des- 
cribed in  Section  3.2.  In  that  section  it  \vas  determined  that  the  best  choice 
of  q was  related  to  the  rms  value  of  the  acceleration  rate  for  a particular 
trajectory.  For  a highly  evasive  trajectory,  the  rms  value  of  a^  is  large  and 
a large  value  of  q is  required  to  minimize  the  rms  estimation  error.  This 
effectively  provides  a large  filter  bandwidth  permitting  the  filter  estimates  to 
track  the  changing  target  acceleration.  LikeAvise  for  a nonevasive  trajectory 
where  the  acceleration  rate  is  small,  q should  be  small  so  that  the  filter  has  a 
low  bandwidth  to  achieve  good  noise  suppression.  In  the  adaptive  design  of 
Fig.  3.3-1,  q is  calculated  on-line  using  the  estimates  of  a^  from  the  fourth- 
order  filter,  in  the  manner  described  below. 

Figures  3.3-2  and  3.3-3  show  the  results  of  a series  of  simulations 

made  to  find  the  optimum  choice  of  the  design  parameter  assuming  the  target 

acceleration  rate  is  known.  Figure  3.3-2  is  a plot  of  the  steady  state  mean 

square  normal  acceleration  estimation  error  produced  by  a third-order  fixed 

filter  which  is  tracking  a target  whose  acceleration  rate  is  held  constant.  Each 

curve  is  for  one  particular  value  of  a and  shows  how  the  estimation  error 

0 

varies  with  the  choice  of  q.  Except  for  the  fact  that  acceleration  rate  is  con- 
stant, the  simulation  conditions  here  are  the  same  as  those  specified  in  Section 

3.1,  Eq.  (3.1-7).  It  is  clear  from  Fig.  3.3-2  that  for  each  value  of  a^ 

6 
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Fip:ure  3.3-2 


Mean  Snuare  Normal  Acceleration  Estimation  Error 
versus  q k r Various  Constant  Values  of  Acceleration 
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APPROXIMATE 


Figure  3.3-3 


NORMAL  ACCELERATION  RATE,  a, 

Optimum  Design  Parameter  versus  Target 
Acceleration  Rate:  Derived  from  Figo  3.3-2 


there  is  an  optimum  choice  of  q.  The  dashed  Ime  in  the  figure  is  the  locus  of 
the  minima  of  the  curves.  The  solid  curve  in  Fig.  3. 3-3  is  a plot  of  the  best 
value  of  the  design  parameter,  designated  q*,  for  each  value  of  acceleration 
rate.  This  shows  how  to  choose  q for  the  variable  bandwidth  filter  if  the 
target  acceleration  rate  is  known.  To  simplify  the  implementation,  the  dashed 
curve  approximation  in  Fig.  3.3-3  is  used,  * ' ’ch  is  described  by  the  equation 


q*  = 4,8  {k/ 


(3.3-1) 


This  approximation  is  justified  because  it  results  in  only  small  errors  in  the 
choice  of  q*  ; much  larger  errors  will  be  incurred  in  estimating  the  value  of 
a^  to  be  used  in  Eq.  (3.3-1).  The  fact  that  q*  depends  upon  the  square  of  the 
acceleration  rate  is  in  agreement  with  the  general  relationship  between  q and 
the  rms  value  of  a indicated  in  the  experiments  of  Section  3.2. 

u 
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Replacing  a.  by  a.  in  Eq„  (3.3-1),  the  third-order  filter  bandwidth 

P P 

parameter  is  adjusted  according  to 

q = 4.8  (3.3-2) 

and  the  filter  gains  k^,  needed  to  mechanize  Eq.  (3.2-4),  are  computed  on-line 
as  shown  in  Fig.  3.3-4. 

R.11915 

THIRD-ORDER 

ADAPTIVE  MATRIX  RICCATI  VARIABLE  BANDWIDTH 

ALGORITHM  EQUATION  FILTER 


Figure  3.3-4  Functional  Diagram  of  Adaptive  Bandwidth  Filter  Design 

The  design  in  Fig.  3.3-4  is  based  upon  the  specific  set  of  e:q)eri- 
mental  conditions  defined  in  Eq.  (3.1-7).  Perhaps  the  quantity  that  is  least 
accurately  known  is  the  rms  value  of  the  measurement  noise.  If  a different 
level  of  noiSv^^  were  used  in  the  reference  model  then  Eq.  (3.3-2)  would  have  to 
be  rederived.  Thus,  this  adaptive  design  will  be  sensitive  to  inaccurate  know- 
ledge of  the  measurement  noise  statistics. 

Finally,  the  bandwidth  parameter,  q,  for  the  fixed  fourth-o;  der  f ‘Iter 
must  be  selected.  Si»-?e  this  quantity  is  to  remain  fixed  and  its  optimum  value 
depends  on  the  rms  acceleration  rate,  there  is  no  a priori  ’’best”  choice. 

Here  its  value  is  arbitrarily  picked  to  yield  the  best  estimates  of  a^  for  the 
nominal  trajectory  in  Fig.  2.3-2(b).  This  value  is  deterinined  from  Fig.  3.2-6 
to  be  approximately  7.5  ft  /sec  . Tliis  completes  the  choice  of  design  par- 
ameters for  the  adaptive  bandwidth  filter  to  be  used  in  the  simulations  described 


in  the  next  cliapter,  where  the  various  adaptive  techniques  are  compared. 
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3.3.2  Hypothesis  Testing;  Filter 

This  section  presents  an  application  of  the  general  hypothesis 
testing  theory  discussed  in  Appendix  B.  The  hypothesis  testing  filter  partially 
circumvents  the  problems  of  fixed  filter  designs  by  optimizing  the  filter  design 
parameters  on-line.  It  assumes  ^*hat  over  some  time  interval,  T,  the  optimum 
filter  is  one  member  of  a set  of  N possible  known  filters  — i.e. , is  one  of  N 
hypotheses.  These  filters  are  run  in  parallel  and  an  algorithm  operates  on 
each  set  of  estimates  to  calculate  the  probability  that  each  hypothesis  is  cor- 
rect. The  final  state  estimate  is  the  sum  of  the  estimates  from  all  N filters, 
each  weighted  by  the  probability  that  the  corresponding  filter  is  optimum. 

Figure  3.3-5  shows  the  structure  of  the  hypothesis  testing  filter  for 
the  target  tracking  problem.  For  the  purpose  of  limiting  the  amount  of  com- 
putation required,  only  three  parallel  filters  are  considered,  each  based  on 
the  third-order  model  developed  in  Section  3.2.  Different  values  of  the  design 


FILTER  1 R-11642 


Figure  3.3-5  Hypothesis  Testing  Filler 
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parameter  q are  chosen  to  represent  different  hypotheses  about  the  target 
behavior.  The  low,  medium,  and  high  values  of  q shown  in  the  figure  cor- 
respond to  low,  moderate,  and  high  target  evasive  xnaneuver  levels. 


Followuig  the  notation  of  Appendix  B,  H is  the  hypothesis  that  Filter 

2 3 

1 is  optimum  over  the  interval  T.  H and  H are  designated  in  a similar 

manner.  The  interval  measurement  history,  Z^,  up  to  and  including  the  most 

recent  measurement  z is  the  set 

®n 


■ V ^ •••  ’ \ 

I n-m-l 


n-l  \ 


where  m is  the  number  of  measurements  taken  in  the  interval. 


Given  the  three  estimates  of  the  system  state  x during  the  interval, 
denoted  by  x^,  x^ , and  x and  generated  by  Kalman  filters  designed  accord- 
ing to  the  procedure  described  in  Section  3. 2,  the  probability  of  each  hypothesis 
conditioned  on  the  measurements,  denoted  by  p {H^|Z  ( , is  given  recursively 
by 


pIhV„I  = 


n-l  M 


S Vi)i>{«''lVi} 

k = 1 


(3.3-3) 


where  f | is  a knouai  normal  probability  density  function  defined 

in  Appendix  B.  In  this  manner,  the  probability  that  each  hypothesis  is  true  is 
calculated  at  every  sampling  time  in  the  interval. 

The  estimate  produced  by  (he  hypothesis  testing  filter  immediately 
after  a measurement  is 
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3 

ia<*l  • S 1 (3.3-4) 

k=  1 

with  the  corresponding  covariance  matrix  given  by 

3 

■’nM  ■ S (3.3-5) 

k=  1 


To  complete  the  filter  design,  it  is  necessary  to  select  the  reset 
interval,  T,  and  the  a priori  probabilities  needed  to  initialize  Eq.  (3.3-3)  at 
the  beginning  of  each  intervaL  The  interval  T should  reflect  tlie  approximate 
length  of  time  over  which  the  target’s  behavior  is  likely  to  remain  constant. 

The  a priori  probabilities  for  each  hypothesis  should  be  chosen  using  any 
knowledge  available  which  indicates  the  relative  probabilities  of  the  various 
types  of  target  behavior.  For  example,  if  it  is  reasonably  certain  that  the 
target  will  have  a highly  evasive  trajectory,  then  the  a priori  probability  of 
the  third  h>q>othesis  should  be  set  much  higher  than  that  of  the  others. 

At  the  end  of  each  interval  the  filters  are  reset  to  test  for  the  possi- 
bility that  a change  has  occured  in  the  target  manuver  characteristics.  This 
is  done  by  resetting  all  tlircc  estimates  and  covariance  matrices  to  the  com- 
bined estimate  and  combined  covariance  matri.x  i s calculated  by  Eqs.  (3.3-4) 
and  (3.3-5)  at  the  end  of  the  interval,  and  the  probability  of  each  hypothesis  is 
reset  to  its  a priori  value.  This  procedure  allows  the  tracking  filter  to  periodi- 
cally adapt  to  new  target  behavior. 

Tliis  design  has  i\vo  disadvantages.  First,  it  requires  approximately 
three  times  as  much  compulation  as  the  third-order  fixed  filler  discussed  in 
Section  3.2  and  about  30Vo  more  compuiation  than  the  adaptive  bar.dwidth  design 


3-2G 


L 


L 

L 


( 


! 

i 

! 

I 


THE  ANALYTIC  SCIENCES  CORPORATION 


which  only  requires  that  two  filters  be  implemented.  Second,  it  will  be  demon- 
strated in  the  results  of  Chapter  4 that  this  desip?i  is  very  sensitive  to  the 
proper  choice  of  the  design  value  for  the  rms  measurement  noise  level;  that  is, 

If  the  noise  level  is  not  knouTi  accurately  and  the  filter  is  designed  assuming 
the  wrong  value,  the  estimation  accuracy  is  seriously  degraded.  This  limits 
the  use  of  the  hypothesis  testing  fiUer  to  tracking  situations  where  the  meas- 
urement noise  level  is  fairly  well  known  ? priori.  Otherwise,  this  design  is 
theoretically  sound  and  can  be  expected  to  u^rk  quite  well  as  long  as  the  assump- 
tions upon  which  it  is  based  are  not  seriously  violated. 

3.3.3  Residual  Testing  Filter 

1 

j 

The  residual  testing  filter  is  similar  to  the  design  discussed  in  the  | 

last  section;  its  complete  details  arc  given  in  Appendix  B.  The  motivation  for  i 

this  particular  design  is  the  noise  sensitivity  problem  discovered  when  the  j 

hypothesis  testing  filter  was  tested  for  various  rms  measurement  noise  levels;  j 

it  is  an  attempt  to  retain  the  advantages  of  the  hypothesis  testing  scheme  while 
reducing  its  noise  sensitivity.  Figure  3.3-6  shows  the  structure  of  the  residual 
testing  filter,  where  the  individual  filters  1,  2,  and  3 are  exactly  the  same 
design  as  those  in  Fig.  3.3-5.  That  is,  Fillers  1,  2,  and  3 are  designed  for 
low,  moderate,  and  high  target  maneuver  levels,  respectively.  The  difference 
between  this  approach  and  the  hypothesis  testing  filler  is  the  manner  in  which 
the  final  slate  estimate  is  calculated  from  the  three  individual  stale  estimates. 

This  residual  testing  tecluiique  docs  not  attempt  to  calculatve  the 
probability  that  one  of  the  individual  filters  is  optimum  over  an  interval. 

Instead  it  observes  each  state  estimate  over  the  interval  T,  at  the  end  of  which 
a decision  is  made  as  to  which  estimate  is  best.  Over  the  next  interval,  the 
residual  testing  filter's  output  is  the  output  of  that  individual  filter  which  worked 
best  over  the  previous  interval. 
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Figure  3.3-6  Residual  Testing  Filter 


The  best  of  the  three  parallel  filters  over  an  interval  is  defined  to  be 
that  one  whose  residual  sequence  has  the  smallest  mean  square  value.  Elements 
of  the  residual  sequence  for  the  filter  are  defined  as 

'i;  ■ \ - "ii:  <-> 

The  mean  square  value  is  computed  from  the  residual  sequence  using 


n +M-1 
o 


n = n 

u 


where  n^  is  the  \’alue  of  the  time  indox  n at  the  beginning  of  the  interval,  T, 
and  M is  the  total  number  of  samples  in  the  interval.  This  calculation  is  per- 
formed for  each  filter,  i.c.,  k - 1,  2,  and  3,  and  at  iUo  end  of  the  itdcrval  thai 
filter  which  has  the  lowest  value  of  g*^  is  judged  best.  Its  estimates  are  the 
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output  of  the  residual  testing:  filter  over  the  next  interval.  The  justification  for 
this  adaptive  procedure  is  discussed  in  detail  in  Appendix  B. 

At  the  end  of  every  interval,  each  filter  is  reset  to  the  state  estimate 
and  covariance  matrix  of  the  filter  which  is  judged  to  be  best.  This  is  done  to 
permit  adaptation  to  changes  in  target  behavior.  The  design  of  the  residual 
testing  filter  is  completed  by  choosing  the  reset  interval,  T.  This  must  be  a 
practical  compromise  between  the  objective  of  achieving  rapid  adaptation, 
which  requires  a small  T,  and  achieving  an  accurate  estimate  of  the  rms  level 
of  the  residuals,  which  requires  a reasonable  averaging  interval  in  Eq.  (3.3-7). 

The  computational  requirements  of  the  residual  testing  scheme  are 
essentially  the  same  as  those  for  the  hypothesis  testing  filter.  However,  it  has 
the  advantage  that  it  is  less  sensitive  to  differences  between  the  actual  and 
design  values  for  the  rms  measurement  noise  level  than  the  other  adaptive 
filter  design  techniques  as  demonstrated  by  the  simulation  results  presented 
m the  next  chapter. 

In  this  chapter,  practical-suboptimum  target  tracking  filters  were 
designed  and  discussed  in  relation  to  the  optimum  design.  Both  fixed  and 
adaptive  filter  structures  applicable  to  the  target  tracking  problem  were  con- 
sidered. In  the  next  chapter,  these  designs  are  compared  and  evaluated  in  a 
variety  of  tracking  situations  to  see  which  of  them  would  make  good  target 
tracking  fillers  in  a practical  application. 


This  cliapter  gives  the  results  of  simulations  designed  to  compare 
the  performance  of  the  fixed  and  adaptive  tracking  filters  discussed  in  the  pre- 
ceeding  chapter.  In  the  first  half  of  this  chapter  the  results  are  presented  in 
terms  of  the  accuracy  with  which  the  filter  can  estimate  the  target’s  position, 
velocity,  and  acceleration  and  in  the  second  half  these  same  results  are  pre- 
sented in  terms  of  the  accuracy  with  which  the  filter  can  predict  the  future  posi 
tion  of  the  target.  The  filters  are  evaluated  for  different  levels  of  target 
evasive  behavior,  levels  of  measurement  noise,  data  rates,  amounts  of  correla 
tion  in  the  measurement  noise,  and  target  ranges.  In  Section  4.4,  tho  sensi- 
tivity of  prediction  error  to  projectile  velocity  and  target  closing  velocity  is 
also  discussed.  The  advantages  and  disadvantages  of  each  filter  are  analyzed 
and  related  to  the  specific  properties  of  each  design. 


4. 1 TRACKING  FILTER  EVALUATION  PROCEDURE 


To  give  a fair  and  logical  comparison  between  several  filters,  it  is 
necessary  to  define  a nominal  tracking  situation.  Each  algorithm  is  designed 
based  upon  this  nominal  situation;  performance  results  are  then  compared  for 
a variety  of  different  operating  conditions.  The  following  set  of  conditions 
define  the  nominal  tracking  situation: 


The  target’s  normal  acceleration  profile  is  *;iven  in 
Fig.  2.3-2(b),  the  nominal  case  with  the  rms  value 
of  acceleration  rate  (a  ) equal  to  30  ft/sec^. 
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L. 


The  measurement  noise  sequence  v (t  ) is  a white 

u n 

gaussian  sequence  with  an  rms  level  of  = 1.4 

mrad,  i.e.,  the  rms  angular  measurement  error 
is  1.4  mrad,  or  equivalently,  the  rms  position 
measurement  error  normal  to  the  line-of-sight  at 
10, 000  ft  is  14  ft. 

The  data  rate  is  10  measurements  per  second. 

The  initial  estimate  at  t = 0 for  all  filters  is  = 0. 

The  initial  covariance  of  the  estimation  error  at 
t = 0 for  all  filters  is 


"o'-’ 


2.0xl0“®(rad)^ 

0 

0 


0 0 
7.lxl0“^(rad/8ec)*  0 

0 2,7xlO^(fl/»ccV 


• The  initi.vl  condition  on  the  target  is 


) 

'2.0  X 10"^ 

(rad) 

’e(o)  ■ 

x(0)  = 

I 

4.0x10"^ 

(rad/sec) 

- 

e(0) 

I. 

.2.0x10^ 

(ft/sec^)  _ 

This  nominal  tracking  situation  is  selected  to  be  generally  representative  of 
digital  gun  fire  control  system  capabilities.  The  tracking  filters  are  evaluated 
by  holding  all  of  the  above  conditions  constant  except  one,  which  is  varied  to 
test  the  sensitivity  of  each  design  to  a change  in  that  condition  only.  Six  dif- 
ferent tracking  filter  designs  were  investigated.  These  are  specified  below. 

Filter  A --  Third-Order  Fixed  - Filter  A is  the  third-order  fixed 
design  discussed  in  Section  3.2  which  models  the  larget’s  normal  acceleration 
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as  the  output  of  a single  integrator  driven  by  white  gaussian  noise  with  a 
spectral  density  of  q.  For  this  design,  q is  chosen  to  be  equal  to  750  (ft  /sec  ), 
which  is  approximately  the  best  value  for  the  nominal  target  acceleration  pro- 
file, as  indicated  by  Fig.  3.2-1. 


Filter  B — Adaptive  Bandwidth  - The  adaptive  bandwidth  filter  is  dis- 
cussed in  Section  3.3.1.  Its  structure  is  specified  in  Figs.  3.3-i  and  3.3-4. 

Filters  CE,  CL,  and  CH  — Hypothesis  Testing  - The  design  of  the 
hypothesis  testing  filter  is  discussed  in  Section  3.3.2  and  its  structure  is 
shown  in  Fig.  3.3-5.  In  each  case  the  reset  interval  chosen  is  T = 2 seconds. 
This  choice  is  made  because  any  shorter  interval  does  not  give  the  probability 
calculation  algorithm,  (Eq.  (3.3-3)),  sufficient  time  to  converge  to  the  correct 
hypothesis  in  the  nominal  tracking  case,  whereas  a longer  reset  time  reduces 
the  filter’s  ability  to  adapt  rapidly  if  the  tracking  situation  changes.  Filters 
CE,  CL,  and  CH  are  exactly  the  same  design  except  for  the  choice  of  the  a 
priori  probabilities,  p )h^{  , at  the  beginning  of  each  reset  interval.  Filter 
CE  assumes  that  all  three  of  the  hypotheses  are  equally  probable  at  tha  start 
of  each  interval.  Filter  CL  assumes  that  there  is  a high  probability  that  the 
low  maneuver  filter  is  optimum;  specifically 


p{hM  = 0.90 
pjH^j  = pjH^I  = 0.05 


Filer  CL 


Filter  CH  assumes  that  there  is  a high  probability  that  the  filter  designed  for 
the  highly  evasive  manuver  case  is  optimum.  The  a priori  probabilities  chosen 
are 


P(kS  = p(H^)  = 0.05 


f\  e\r\ 


Filter  CH 
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Filter  D Residual  Tcstinp^  - The  design  of  the  residual  testing, 
filter  is  discussed  in  Section  3.3.3  and  its  structure  is  shown  in  Fig.  3.3-6. 
Here  again  the  reset  interval  chosen  is  2 seconds  because  this  is  the  shortest 
interval  found  to  permit  a realistic  choice  between  the  three  parallel  filters 
based  on  the  analysis  of  the  residual  sequence.  For  shorter  intervals  there  is 
not  enough  data  to  make  an  accurate  estimate  as  to  which  filter  is  operating 
best  under  the  nominal  tracking  conditions. 


follows: 


The  six  filters  and  their  letter  designations  are  summarized  as 

• A - Third-Order  Fixed  Filter 

• B - Adaptive  Bandwidth  Filter 

• CE  - Hypothesis  Testing  Filter  (Equal  Probability) 

• CL  - Hypothesis  Testing  Filter  (Low  Evasion  Probability) 

• CH  - Hypothesis  Testing  Filter  (High  Evasion  Probability) 

• D - Residual  Testing  Filter 


4. 2 FILTER  PERFORMANCE  AND  SENSITIVITY  COMPARISONS 

In  this  section  the  six  filters  described  above  are  compared  in  terms 
of  rms  errors  in  the  estimates  of  the  current  (not  predicted)  target  position, 
velocity,  and  acceleration  under  different  tracking  situations.  The  rms  esti- 
mation error  shown  on  the  figures  is  calculated  by  taking  the  square  root  of  the 
time  average  of  the  estimation  error  squared  over  the  period  from  twn  seconds 
until  ten  seconds.  This  is  done  to  omit  the  initial  transient  in  the  estimation 
error.  The  transient  is  due  to  the  large  initial  errors  in  the  estimate  s which 
are  rapidly  eliminated  by  all  of  the  filter  designs—  for  example,  see  Fig. 
3.2-3.  In  all  of  the  results  presented,  except  foj  the  range  sensitivity  study, 
only  the  rniS  normal  acceleration  estimation  error  is  showm  because  the 


4-4 


THE  ANALYTJC  SCIENCES  CORPORATION 


accuracies  of  the  angular  position  and  angular  velocity  estimates  for  the 
various  filters  usually  compare  in  the  same  manner.  As  sho\m  in  Section 
4.3,  target  acceleration  estimation  error  is  the  primary  contributor  to  pre- 
dition  error. 


4.2.1  Sensitivity  to  Target  Maneuver  Level 

Figure  4. 2-1  shows  how  each  filter  performs  for  the  three  target 
acceleration  profiles  defined  in  Fig.  2.3-2.  The  results  for  the  six  filters 
are  presented  in  two  parts  for  clarity;  the  results  for  Filters  A and  CE  are 
shown  in  each  part  of  the  figure  for  comparison  purposes.  Note  that  the  rms 
normal  acceleration  estimation  error  goes  up  almost  linearly  as  a function  of 
the  target’s  rms  normal  acceleration  rate  for  each  design..  Thus,  as  the 
target  increases  its  level  of  evasive  maneuvers,  all  designs  have  decreasing 
accuracy. 

For  the  nominal  trajectory— i.e. , the  mildly  evasive  case— the  simple 
third-order  fixed  filter  is  best,  as  expected,  since  it  is  optimized  for  this 
case.  The  other  filters  give  only  a slightly  higher  rms  error,  except  for  the 
hypothesis  testing  filter  (CH)  which  assumes  that  there  is  a high  probability 
that  the  target  will  maneuver  in  an  evasive  manner;  for  the  nominal  trajectory 
this  is  a poor  assumption,  which  is  reflected  in  the  filter's  performance. 

For  the  nonovasive  trajectory,  Filter  CH  again  works  poorly,  as 
expected.  However,  Filters  B and  CL  offer  a signficant  improvement  over 
Filter  A,  the  fixed  third-order  filter.  Filter  B (adaptive  bandwidth)  is  clearly 
best  able  to  recognize  the  fact  that  the  target  is  not  changing  its  normal  ac- 
celeration and  lowers  its  bandwidth  accordingly  by  reducing  q to  achieve  a 
lower  rms  estimation  eri'or.  Filter  CL  also  works  well  because  it  assumes 
that  there  is  a high  probability  that  the  target  is  nonevasive  and,  for  this 
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trajectory,  this  is  an  accurate  assumption*  Filter  CE,  the  equiprobability 
hypothesis  testing  filter,  works  fairly  well  compared  to  the  others  consider- 
ing the  fact  that  it  assumes  that  all  types  of  trajectories  are  equally  probable, 
when  the  actual  trajectory  is  nonevasive.  Filter  D,  the  residual  testing  filter, 
has  only  a slightly  lower  rms  estimation  error  than  that  of  the  much  simpler 
Filter  A*  Figure  4.2-2  shows  time- histories  of  the  acceleration  estimates 
of  the  CE  and  CL  filters  for  the  nonevasive  case  and  demonstrates  the 
importance  of  a priori  probabilities  in  the  hypothesis  testing  design,  since 
the  two  filters  are  exactly  the  same  except  for  the  choice  of  their  a priori 
probabilities. 

For  the  evasive  trajectory,  all  of  the  adaptive  designs  offer  an  im- 
provement over  the  fixed  design.  In  Fig.  4.2-1,  notice  that  the  CH  filter 
gives  the  curve  with  the  smallest  slope.  This  indicates  that  for  sufficiently 
violent  target  manuvers  (more  violent  than  those  tested  here),  the  CH  filter 
would  be  the  best.  Figure  4. 2-3  shows  the  acceleration  estimate  of  the  CE 
filter  for  this  case. 

In  general,  Fig.  4.2-1  sho\vs  that  adaptive  designs  yield  come  im- 
provement over  the  fixed  design  if  the  change  in  target  behavior  is  radical. 

It  also  shows  that  even  the  best  design  has  a high  rms  estimation  error  if  the 
target  behaves  evasively.  Filter  B,  using  the  adaptive  bandwidth  principle, 
gives  the  best  overall  performance.  Tlie  biggest  percentage  improvement  in 
estimation  error  achieved  by  the  adaptive  filters  over  the  fixed  design  is  at 
low  rms  acceleration  rates.  At  high  rates  no  design  significantly  reduces  the 
tracking  error  below  that  of  the  third-order  fixed  filler  which  is  optimized  for 
the  higli  rate  case. 
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(a)  Hypothesis  Testing  Filter  (Equal  Probability)  ~ Filter  CE 


(b)  Hypothesis  Testing  Filter  (Low  Evasion  Proliability)  - Filter  CL 


Figure  4.2-2 


Normal  Acceleration  and  Us  Estimate  using  Hypothesis 
Testing  Filters  CFI  and  CL:  Nonevasivc  Trajectory 
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h“  RMS  ERROR  = 73  ft/sec2  -y 


(OVER  10  see) 


TIME  (sec) 


Figure  4.2-3 


Normal  Acceleration  and  Its  Estimate  Using 
Hypothesis  Testing  Filter  CE:  Highly  Evasive 
Trajectory 


A further  demonstration  of  the  adaptive  designs  at  low  maneuver  levels 
is  provided  in  Fig.  4*.  2-4  where  Filters  B and  CL  are  compared  with  two  ver- 
sions of  the  third-order  fixed  niter  designated  A*  and  A*^.  The  former  is 
desigiK'd  (q  - 7.5  fl^/scc^)  so  thui,  for  the  nonevusive  trajectory,  its  rms 
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mis  NORMAL  ACCELERATION  RATE 


Figure  4. 2-4  Comparison  of  Adaptive  Filters  versus  Fixed 
Filters  Optimi^d  for  the  Nonevasive 
Trajectory 

acceleration  estimation  error  is  one-half  that  achieved  with  Filter  CL.  Thus, 
the  A*  filter  is  very  close  to  the  optimtim  design  (q  = 0)  for  this  trajectory, 
but  its  performance  is  degraded  for  the  nominal  and  highly  evasive  trajectories. 
This  demonstrates  dramatically  the  price  paid  for  a near  optimum  design  if  the 
target  does  not  behave  as  assumed. 

The  filter  designated  A**  is  the  same  as  A*  except  its  value  of  q is 
selected  to  have  the  same  rms  error  as  filter  CL  for  the  nonevasive  case. 

This  fixed  design  does  nearly  as  well  as  the  adaptive  designs  (CL,  B)  for  the 
nominal  trajectory  but  its  performance  is  poor  for  the  highly  evasive  case. 

Thus,  if  the  filters  are  "matched"  to  have  the  same  rms  error  for  a nonevasive 
target,  the  adaftive  filters  perform  much  better  if  the  target  acutally  is  evasive. 
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4.2.2  Sensitivity  to  Measurement  Noise  Level 

All  of  the  filters  have  been  designed  assuming  a nominal  rms  measure-* 
ment  noise  level  of  1.4  mrad.  This  section  describes  what  happens  to  the  rms 
normal  acceleration  estimates  if  the  actual  measurement  noise  level  is  signifi- 
cantly hi^er  or  lower  than  the  design  value.  This  is  an  important  sensitivity 
study  because  the  measurement  errors  in  an  actual  tracking  situation  are  to 
some  extent  determined  by  variable  factors  such  as  weather  and  target  cross- 
section;  thus  the  rms  error  level  may  not  be  known  accurately  at  any  given  time. 
Figure  4. 2-5  shows  the  simulation  results  obtained  for  the  nominal  target 
trajectory. 

For  the  nominal  noise  level,  the  results  are  the  same  as  these  given 
in  Section  4. 2-1  for  the  nominal  trajectory  case.  When  the  noise  level  is 
decreased,  all  the  filters  yield  approximately  the  same  rms  error  except  for 
Filter  D which  gives  the  best  performance.  The  high  noise  case  gives  the 
most  interesting  and  important  sensitivity  results.  Compared  with  Filter  A, 
all  of  the  adaptive  designs  except  Filter  D are  extremely  sensitive  to  meas- 
urement noise  which  is  higher  than  the  assumed  level.  This  sensitivity  to 
measurement  noise  level  is  common  to  many  adaptive  techniques  auvi  motivates 
the  residual  (:esting  scheme,  Filter  D.  This  behavior  will  be  explained  presently. 

The  adaptive  bandwidth  filter  (B)  is  sensitive  to  measurement  noise 
level  because  its  bandwidth  is  controlled  by  estimates  of  acceleration  rate. 

A 

As  the  estimate  a increases,  the  bandwidth  of  the  filter  also  increases 
through  the  mcclianization  of  Eq.  (3.3-2).  This  is  the  desired  result  if  the 
actual  rms  acTceleration  rate  also  increases.  However,  if  the  measurement 
nois^i  level  increases,  then  the  state  estimates  produced  by  the  fourth-order 

A 

fixed  filter  in  Fig.  3.3-1  become  more  noisy;  hence  the  rms  level  of  a.  will 

0 


Figure  4.2-5 


Tracking  Filter  Sensitivity  to 
Measurement  Noise  Level 


4.12 


THE  ANALYTIC  SCIENCES  CORPORATION 


be  higher  due  to  its  increased  noise  content.  Consequently,  the  adaptive  band- 
width filter  cannot  tell  whether  the  estimated  rms  acceleration  rate  has  in- 
creased due  to  an  evasive  maneuver  or  to  higher  measurement  noise,  but  its 
bandwidth  will  be  raised  in  both  cases.  In  the  former  case  its  bandwidth 
should  be  raised;  in  the  latter  it  should  be  lowered.  Thus  the  filter  adapts  in 
the  wrong  direction  when  changes  in  rms  measurement  noise  levels  occur. 


The  hypothesis  testing  filter  suffers  from  the  same  difficulty  des- 
cribed above.  The  design  based  on  the  nominal  noise  level  tells  the  filter  how 
ttuch  noise  to  e^qiect  in  the  estimates.  If  the  actual  noise  level  is  greater  than 
its  assumed  value,  the  filter  effectively  interprets  the  increased  fluctuations  in 
the  measurements  as  caused  by  increased  target  evasive  action  and  errone- 
ously determines  the  high  bandwidth  filter  in  Fig.  3.3-5  to  be  »e  most  probable. 

The  residual  testing  filter  (D)  overcomes  the  noise  sensitivity  problem 
because  its  bandwidth  is  directly  related  to  the  rms  value  of  the  residual 
sequence  of  each  parallel  filter  in  Fig.  3.3-6.  As  discussed  in  Section  3.3-3 
and  Appendix  B.  2,  this  sequence  is  a direct  reflection  of  the  actual  estimation 
error.  Filter  D determines  which  of  the  parallel  filters  has  the  lowest  rms 
residuals  and  therefore  provides  the  smallest  estimation  error.  Another  way  of 
looking  at  this  is  to  notice  that  the  hypothesis  testing  and  adaptive  bandwidth 
filters  have  designs  in  which  the  adaptive  algorithm  depends  upon  the  assumed 
level  of  noise.  If  this  level  is  incorrect,  they  adapt  in  the  wrong  direction. 

On  the  other  hand,  the  adaptive  part  of  the  residual  testing  filter  is  independent 
of  the  actual  measurement  noise  level;  therefore  it  adapts  correctly. 

One  possible  means  of  reducing  the  noise  sensitivity  of  Filters  B,  CH, 
CL,  and  CE  would  be  to  design  them  assuming  the  highest  level  of  measurement 
noise  that  might  ever  be  encountered.  This  would  work  well  in  the  high  noise 
case,  but  the  estimates  would  compare  poorly  with  those  of  Filters  A and  D if 


i 

i 
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tb  ^xtual  noise  level  were  much  lower.  In  situations  where  the  measurement 
noise  level  is  poorly  known  or  changing,  a better  approach  is  to  use  a filter 
which  is  designed  for  the  most  likely  noise  level  but  is  insensitive  to  changes 
around  this  design  value;  for  example,  the  residual  testing  filter  designed  for 
the  nominal  noise  level. 


4.2.3  Sensitivity  to  Correlated  Measurement  Noise 

In  the  tracking  problem,  the  sensor  (c.g. , radar)  receiver  noise  is 
usually  wide-band  so  that  white  gaussian  noise  is  a good  model.  However  the 
total  measurement  error  can  also  contain  low  frequency,  correlated  noise 
components  caused  by  a slowly  varying  target  cross  section,  due  to  the  target*s 
own  motion.  The  simulation  described  here  ij  conducted  to  see  how  much  cor* 
related  noise  affects  estimation  accuracy.  Except  for  the  measurement  noise, 
the  tracking  situation  is  the  nominal  case  defined  in  Section  4.1.  Only  Filter 
A is  tested,  but  similar  results  would  be  obtained  for  all  of  the  filters. 

The  measurement  noise  sequence  is  generated  by  driving  the  difference 
equation  for  a low-pass  filter  with  the  output  of  a gaussian  random  number 
generator.  The  details  of  the  measurement  noise  model  used  in  this  simulation 
are  presented  in  Section  2.3,  where  the  correlation  time  is  denoted  by  the 
var  ictu  .or. 

Figure  4.2-6  shows  that  increasing  the  correlation  initially  in- 
creases the  error.  As  the  correlation  time  becomes  very  large,  (i.e.,  t>1 
sec.),  the  error  begins  to  decrease  because  very  low-frequency  noise  has 
little  effect  on  the  estimation  of  a time-derivative  such  as  target  acceleration. 
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CORRELATION  TIME,  T(sec) 


Figure  4. 2-6 


Sensitivity  of  Filter  A to  Measurement 
Noise  Correlation 


4. 2. 4 Sensitivity  to  Data  Rate 


One  means  of  improving  estimation  accuracy  is  to  increase  the  infor- 
mation available  to  the  filter  by  increasing  the  frequency  with  which  measure- 
ments are  taken.  Thus  far  ten  measurements  per  second  has  been  assumed  to 
be  a reasonable  data  rate  for  the  tracking  system.  To  see  how  much  improve- 
ment is  possible,  Filter  A was  exercised  at  data  rates  of  twenty  and  fourty 
measurements  per  second;  the  results  are  shown  in  Fig.  4. 2-7.  It  is  clear 
from  this  gnq>h  that  a hi^r  data  rate  will  improve  estimation  accuracy; 
However,  this  improvement  is  achieved  at  the  ejqpense  of  increased  computer 
time  requirements— i.e.,  if  the  data  rate  is  doubled,  then  it  takes  approxi- 
mately twice  as  much  computer  time  to  calculate  the  estimates  for  a given 
computer  design;  thus;  a faster  computer  might  be  required.  Data  compression 
^refiltering)  techniques  could  probably  be  useful  for  reducing  the  computation 
time  required  at  high  data  rates  (Ref.  20). 


ANALYTIC  SCIENCES  CORPOF^ATION 


Figure  4. 2-7  Sensitivity  of  Filter  A to  Data  Rate 


4.2.5  Sensitivity  to  Target  Range 


All  of  l«e  resiiUs  prewnted  la''  1“™  assumed  tl»t  the  target  ts  at  a 
constant  ranee  ol  10, 000  ft.  However,  an  attacking  target  wlU  be  at  this  rMge 
only  once  during  Its  trajectory;  thus  It  Is  Important  to  know  how  the  estimation 

errors  of  ecich  filter  vary  with  range. 


H,  determine  the  estimation  error  sensitivity  to  range,  a set  of 
simulations  was  performed  under  the  nominal  tracking  situation  with  the  range 
parameter  r changed  for  Filters  A,  B,  CE,  mid  D.  Figure  4.2-8  shows  msutts 
tor  ranges  tmm  2,500  to  25,000  ft.  hi  part  (a)  of  this  figure  the  rms  position 
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estimation  error  in  feet  is  plotted  against  range  for  the  four  filters,  This 
error  is  calculated  using  the  small  angle  approximation  for  the  angle  estimation 
error.  Thus  the  position  estimation  error  e^  is  given  by 

Cp  = r sin  (e^)  = re^ 

where  r is  range  in  feet  and  e^  is  angular  estimation  error  in  radians.  Like- 

V 

wise  in  part  (b)  of  the  figure,  the  velocity  estimation  error  in  feet  {>er  second 
is  found  by  using 


where  e^  is  equal  to  0 - the  angular  velocity  estimation  error.  In  part  (c) 
of  Fig,  4.2-8,  the  rms  normal  acceleration  estimation  error  is  shown. 

As  expected,  the  estimates  become  more  accurate  as  the  range  of  the 
target  decreases  since  this  is  an  angle  measuring  system.  For  a constant  rms 
angular  measurement  error,  the  measured  position  error  decreases  linearly 
with  range,  and  consequently  the  estimated  position,  velocity,  and  acceleration 
become  more  accurate  as  the  signal-to-noise  mtio  of  the  data  improves.  In 
practice,  the  rms  angular  measurement  error  may  vary  with  range;  this  effect 
is  not  included  in  this  study  because  it  is  dependent  upon  the  particular  target 
encountered  and  tracking  sensor  used. 


All  of  the  filters  have  approximately  the  same  range  sensitivity  as 
shown  in  Fig.  4. 2-8  except  for  the  hypothesis  testing  filter  which  does  poorly 
at  the  longer  ranges.  This  seems  to  be  due  to  the  fact  that  the  higher  effective 
noise  level  at  longer  ranges  (i.e.,  the  measured  position  error  is  larger) 
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4.2,6  Summary  of  Filter  Estimation  Accuracy 

The  results  described  in  this  section  show  the  merits  of  each  adaptive 
design  in  terms  of  estimation  accuracy.  If  the  rms  measurement  noise  level 
is  well  known,  then  all  of  the  adaptive  designs  offer  some  advantage  over  the 
fixed  filter  (A)  for  targets  with  w^'dely  different  maneuvermg  characteristics. 
For  this  case  the  adaptive  bandwidth  filter  (B)  appears  to  give  the  best  tradeoff 
between  conqplexity  and  estimation  accuracy.  The  price  paid  for  this  design  is 
the  requirement  that  both  a third-order  and  a fourth-order  tracking  filter  be 
implemented,  giving  a more  than  two-fold  increase  in  computational  require- 
ments over  the  third-order  fixed  filter.  If  accurate  a priori  knowledge  is 
available  about  the  probabilities  of  various  target  trajectories,  then  the  hypo- 
thesis testing  filter  (CE,  CL,  or  CH)  is  the  best  choice.  However  it  is  more 
complex  than  Filter  B and  it  is  especially  sensitive  to  changes  in  measurement 
noise  level  and  target  range,  so  it  must  be  used  with  caution. 

For  tracking  situations  where  the  measurement  noise  level  is  unknown 
or  changing,  the  fixed  third-order  filter  and  the  residual  testing  filter  yield 
the  best  performance.  The  residual  testing  filter  can  offer  better  estimates 
than  the  fixed  filter  if  the  target  changes  its  behavior  over  a wide  range,  or  if 
the  measurement  noise  Traries  si:ifficiently.  The  cost  of  this  is  a three-fold 
increase  in  computation  over  the  fixed  third-order  design. 


4. 3 PREDICTION  ERROR  EVALUATION  PROCEDURE 


In  this  section,  the  procedure  used  to  evaluate  the  tracking  filter 
designs  in  term'i  of  prediction  error  is  qsecified.  Prediction  error  is  defined 
as  the  error  made  by- the  filter  in  predicting  the  future  position  of  the  target. 

A simple,  but  realistic,  scenario  for  an  attacking  target  provides  the  basis  for 
plots  of  rms  predicted  position  error  versus  time-to-go,  t^.  The  position 
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error  is  measured  normal  to  the  true  line-of-sight  and  time-to-go  is  defined  as 
the  time  remaining  until  the  attacking  target  reaches  the  ship.  This  is  a more 
practical  basis  for  tracking  filter  evaluation  than  comparing  target  trajectory 
estimation  errors  because  the  ultimate  effectiveness  of  a gun  fire  control 
system  is  determined  by  the  projectile  miss  distance  which  depends  upon  the 
rms  predicted  target  position  error.  Other  sources  of  miss  distance  such  as 
ballistic  dispersion  and  gun  pointing  error  are  not  considered  here*  >^owever, 
these  will  tend  to  be  much  smaller  than  the  target  prediction  errors  for  highly 
evasive  targets. 


The  simplified  scenario  for  an  attacking  target  depicted  in  Fig.  4.3-1 
is  based  upon  several  assumptions^  which  are  not  necessary  in  an  actual  situa- 
tion,but  make  possible  a clear  comparison  of  the  performance  of  different 
tracking  filters.  First,  assume  the  target  is  approaching  the  ship  at  a constant 
radial  or  closing  velocity  denoted  by  v^,  and  its  evasive  action  consists  of 
motion  perpendicular  to  the  line-of-sight.  Thus,  the  target's  range  as  a func- 
tion of  time  is  assumed  known,  but  its  angular  position,  6,  must  be  estimated 
and  predicted.  Second,  it  is  assumed  that  the  projectile  fired  from  the  ship 
tiavels  with  constant  radial  velocity,  Vp.  If  the  target  is  at  range  r when  it  is 
fired,  then  using  the  above  assumptions  the  range  at  intercept,  r^  defined  as 
the  range  at  which  the  projectile  and  target  are  the  same  distance  from  the  ship, 
is  given  by 


= 


Vp  + Vt 


and  the  time-to-go,  t 


go 


is 


(4.3-1) 


(4.3-2) 
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SHIP 

Figure  4. 3-1  Target  Prediction  Scenario 


The  prediction  time,  t^,  is  defined  as  the  time  required  for  the  projectile  and 
the  target  to  reach  the  intercept  range  r^  measured  from  the  time  when  the 
projectile  is  fired.  Thus  t is  given  by 

Mr 


or 


(4.3-3) 


r 


Vp+VT 


(4.3-4) 


Using  Eqs.  (4.3-1)  through  (4.3-4),  range,  intercept  range  and  prediction 
time  as  functions  of  time-to-go  are 


(4.3-5) 
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r = 

I Vp  + V^ 


(4.3-6) 


t = 


Vp  + V^ 


(4.3-7) 


These  three  equations  along  with  choices  for  Vp,  v,p,  and  a range  of  values 
for  tg^  define  the  variables  needed  for  target  prediction.  The  final  step  is  to 
find  the  rms  predicted  position  error  as  a function  of  time-to-go. 

For  a particular  value  of  t^^,  the  target  is  at  range  r and  the  track- 
ing filter  provides  estimates  of  its  position,  velocity  and  acceleration  normal 
to  the  line-of-sight.  These  estimates  are  then  used  to  predict  the  position  of 
the  target  t seconds  later  when  it  reaches  the  intercept  range.  The  difference 
between  the  predicted  and  actual  target  positions  is  the  predicted  position  error 
given  by 


r r 

'OD  • *'i  'VT 


(4.3-8) 


where  e^^  is  the  predicted  position  error  perpendicular  to  the  line-of-sight  at 
the  intercept  range,  e^Ct^),  e^(t^),  e^  (t^),  and  e^  (t^)  are  the  position,  velocity, 

acceleration,  and  acceleration  rate  estimation  errors  at  the  start  of  prediction, 
and  t^  is  the  prediction  time.  Equation  (4.3-8)  follows  directly  from  the  third- 
order  target  model  assumed  in  designing  the  filters  and  the  assumption  that  the 
acceleration  rate  is  constant  during  the  prediction  interval.  Figure  4.3-2  shows 
how  the  last  three  terms  of  Eq.  (4.3-8)  affect  the  error  in  predicted  position  as 
a function  of  t^,  for  various  values  of  the  estimation  errors.  Evidently  the 
predicted  posidon  error  increases  rapidly  with  prediction  time  if  significant 
errors  exist  in  the  estimates  of  target  acceleration  and  acceleration  rate. 
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rRCDiCIION  TIMt,  I («M) 


(a)  Prediction  Error  Due  to  versus  Prediction  Time 


HICOIC1IOI4  1IIK.  liKl 


(b)  Prediction  Error  Due  to  e.  (Iq)  versus  Prediction  Time 

V 


^^*,t 


(c)  Prediction  Error  Due  to  «?•  (;q)  versus  Prediction  Time 

Figure  4. 3-2  Effect  of  Estimation  Error  Components  on  Prediction 
Error  versus  Prediction  Time 
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V . 


Because  the  estimation  errors  in  Eq.  (4.3-8)  are  random,  prediction 
accuracy  can  be  described  in  terms  of  the  rms  value  of  the  predicted  position 


error  as  a function  of  prediction  time,  denoted  by  *uid  defined  by 


i-aj 


(4.3-9) 


where  E [.]  denotes  expected  value.  The  expansion  of  Eq.  (4.3-9)  with  sub- 


stitution from  Eq.  (4.3-8)  leads  to  an  expression  for  which  contains  cross 
terms  of  the  form  E [e  e ],  E [cpe^  ],  etc.  Recall  that  in  Section  3.1,  it  is 
shown  that  the  position,  velocity,  and  acceleration  estimation  errors  are  highly 


correlated.  In  this  case,  a good  approximation  is 


® [VaJ  ■ H'J]  *^  [<])' 


(4.3-10) 


Similar  e:q)ressions  apply  for  the  other  croi.>s-product  terms.  Substituting 
the  latter  into  Eq.  (4.3-9)  produces 


2 3 


(4.3-11) 


where  (t^)  is  the  rms  value  of  eic.  This  approximation  is  worst  case 

in  the  sense  that  less  than  perfect  correlation  between  the  estimation  errors 
leads  to  an  actual  value  of  rms  prediction  error  less  th^n  that  given  by 
Eq.  (4.3-11). 


To  calculate  the  rms  value  of  e using  Eq.  (4.3-11)  it  is  necessarv 

rH 


to  know  the  prediction  time,  t^,  and  the  various  ims  estimation  errors  at  the 


start  of  the  prediction  interval,  which  are  functions  of  the  range  r.  The  quan 


tities  tp  and  r are  given  by  Eqs.  (4.3-5)  and  (4.3-7)  respectively  as  functions 


of  time-to-go,  t . The  range  sensitivity  studies  in  Section  4.2-5  indicate 
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how  the  rms  position,  velocity,  and  acceleration  estimation  errors  vary  with 
range.  The  rms  acceleration  rate  estimation  error,  a*  , is  equal  to  the  actual 

. 

rms  target  acceleration  rate;  i.e. , the  rms  value  of  a. , because  the  filters  do 

V 

not  estimate  a^.  It  was  found  in  Section  3.2  that  attempts  to  estimate  this 
variable  using  a fourth -order  filter  actually  degrade  the  overall  estimation 
accuracy;  the  rate  estimates  are  so  poor  that  they  may  as  well  be  set  equal  to 
. zero.  The  latter  policy  is  followed  here.  In  subsequent  calculations  of  the  rms 
predicted  position  error  this  term  will  be  neglected,  which  is  the  same  as 
assummg  that  the  target  has  a constant  acceleration  over  the  prediction  interval. 
The  effect  of  a change  in  target  acceleration  during  the  interval  can  be  deter- 
mined from  Fig.  4.3-2(c)  and  added  to  the  prediction  error. 

With  the  above  convention  for  handling  a-  , the  rms  predicted  posi- 
tion  error  is  giver*  by 
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(4.3-12) 


The  three  rms  estimation  errors  needed  by  tais  equation  are  calculated  as 
follows:  For  each  filter,  normalized  range  sensitivity  curves  are  plotted  and 
approximated  by  curves  of  the  form  y = shown  in  Fig.  4.3-3,  which  give  a 
good  fit  to  the  normalized  range  sensitivity  curve  of  each  filter.  Figure  4.2-8 
shows  the  original  range  sensitivity  of  position,  velocity,  and  acceleration 
estimation  error  for  Filters  A,  B,  CE,  and  D.  These  curves  were  plotted 
again  with  the  horizontal  scale  normalized  to  10,000  ft— i.e.,  1 corresponds 
to  10,000  ft—  and  the  vertical  scale  normalized  to  the  rms  estimation  error  at 
10,000  ft— i.e. , 1 corresponds  to  a particular  filter’s  rms  position,  velocity, 
or  acceleration  estimation  error  at  10,000  ft.  Figure  4.3-4  is  an  example  of 
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this  for  Filter  A,  showinc;  both  the  approximating;  curves  and  the  actual  simu- 
lation results  taken  from  Fij;.  4.2-8. 

This  procedure  has  been  carried  out  for  each  of  the  tracking  filters 
for  which  range  sensitivity  studies  were  conducted — i.e.,  Filters  A,  B,  CE, 
and  D.  The  approximate  range  sensitivity  curves  obtained  are  based  on  the 
nominal  target  trajectory  and  nominal  rms  measurement  noise  level;  but  they 
can  be  used  to  obtain  the  approximate  rms  estimation  errors  as  a function  of 
range  for  a particular  filter  in  any  tracking  situation  if  the  rms  estimation 
errors  are  known  at  10,  OOG  ft.  Now,  since  time  and  range  are  directly  related 
in  this  scenario  by  the  constant  target  closing  velocity,  the  rms  errors  at  the 
start  of  the  prediction  time  are  defined  by  the  range  at  that  time. 

This  completes  the  detailed  description  of  the  prediction  error  eval- 
uation procedure.  In  summary,  for  a particular  target  closing  velocity  v^ 
and  a particular  projectile  velocity  Vp,  the  prediction  time  t^,  the  target 
range  at  firing  r,  and  the  intercept  range  r^  are  calculated  as  functions 
of  tiine-to-go,  t^.  This  assumes  that  the  target  closes  at  a constant  velocity 
and  evades  only  by  moving  perpendicular  to  the  line-of-sight  and  that  the  pro- 
jectile ti'avels  at  a constant  radial  velocity.  Using  the  rms  estimation  errors 
obtained  for  each  filter  in  each  tracking  situation,  and  the  results  of  the  range 
sensitivity  study  provided  in  Section  4.2,  the  rms  predicted  iX)sition  error 
normal  to  the  line-of-sight  is  calculated  as  a function  of  time-to-go,  assuming 
th:^  target  does  not  change  it.s  acceleration  during  the  prediction  interval.  The 
results  of  the  analysis- -plots  of  rms  predicted  position  error  versus  time-to- 
go  for  each  filter  operating  in  a particular  tracking  situation--arc  presented 
in  Section  4.4. 
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4. 4 PREDICTION  ERROR  PERFORMANCE  AND  SENSITIVITY 

This  section  compares  the  performance  of  Filters  A,  B,  CE,  and  D 
by  sho\dng  how  their  rms  predicted  position  errors  vary  with  time-to-go.  A 
detailed  account  of  how  these  results  are  obtained  is  given  in  Section  4.3 
together  with  an  ejqDlaiiation  of  the  approximations  and  assumptions  involved. 
From  the  results  presented  in  this  se'^^ion  a general  idea  of  the  kill  probability 
could  bo  obtained  by  assuming  a particular  rate  of  fire,  projectile  lethal  radius, 
and  rms  gun  pointing  error.  The  results  not  only  show  the  comparitive  per- 
formance of  the  filters  but  also  indicate  the  overall  accuracy  that  is  likely  to 
be  achieved  with  any  good  tracking  filter  for  the  types  of  target  maneuvers 
considered  in  this  report. 

To  complete  the  specification  of  the  target  attack  scenario  described 
in  the  last  section,  the  projectile  velocity,  Vp,  and  the  target  closing  velocity, 
Vp,  must  be  chosen.  For  most  of  the  results  presented  in  this  section,  these 
quantities  are  taken  i j be 

Vp  = 2600  ft /sec 
Vp  = 1300  ft /sec 

At  the  end  of  this  section,  the  sensitivity  of  the  rms  predicted  position  error  to 
changes  in  these  values  is  shown.  With  this  choice  of  VpOnd  Vp,  the  range  r, 
intercept  range  rj,  and  the  prediction  time  t^  are  completely  specified  as  func- 
tions of  time-to-go  (Eqs.  (4.3-5),  (4.3-C),  and  (4.3-7))  as  shown  in  Fig.  4,4-1. 


4.4.1  Sensitiviix  (o  Target  Maneuver  Level 

The  target  rms  predicted  position  errors  corresponding  to  the  esti- 
mation erurs  acnieved  for  the  three  diffeicnt  target  trajectories,  in  the  otherwise 
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Figure  4.4-1  Target  Range,  Intercept  Range,  and  Prediction  Time 
versus  Time-To-Ck)  for  - 2600  ft  ^sec  and  = 
1300  ft /sec  ^ ^ 


nominal  tracking  situation  defined  in  Section  4.1,  are  given  in  Fig.  4.4-2. 

The  results  are  given  for  Filters  A,  B,  CE,  ami  D,  and  for  ear.:i  of  the  target 
trajectories  defined  in  Fig.  2.3-2.  To  interpiet  these  curves,  assume  that 
only  rins  predicted  position  errors  les.s  than  100  ft  arc  of  interest  and  anythin 
larger  i.s  an  almost  .sure  *'ini.ss**:  llien  the  length  of  time  a particular  curve  is 
belO\V  tins  100  ft  llireshold  is  F4»r  o.'ianmk  . in  Fig.  4.4-2(a)  Cu:  v 


RV3  PRED  CTE^  POSITION  ERROR,  ^ (ft)  RMS  PREOCCTEO  POSITION  ERROR.  (T  (ft)  RMS  PREDICTED  POSITION  ERROR.  <7  ' 
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(a)  Nonevasive  Case  - rms  (a^  ) = 0 
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(b)  Nominal  Case  • rms  {b q)  - 30  ft/scc^ 
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(i;)  Miijiliv  I vusive  Case  nm  (j  j/1  00 


Fij;urc  4A-2  RMS  Prediclod  Position  Krror  versus  Time-To-Go 
Sensitivity  to  Tarj^et  rms  Acceleration  Rate, 

Filters  A,  H,  CL\  and  D 
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B for  the  adaptive  bandwidth  filter  is  below  100  ft  for  approximately  9 seconds 

before  the  target  reaches  the  ship.  Thus  there  are  9 seconds  of  firing  time 

during  which  there  is  some  chance  of  hitting  the  target.  K the  gun  fires  one 

projectile  per  second,  then  there  are  nine  firing  opportunities  before  the  target 

reaches  the  ship.  Suppose  that  one  projectile  is  fired  at  = 9 seconds,  one  is 

fired  at  t =8  seconds  and  so  on  until  t =1  second.  The  rms  predicted 

position  error  associated  with  each  of  these  projectiles  can  be  found  f^Dm  the 

height  of  the  curve  at  the  respective  time  it  is  fired.  Thus  the  rms  error  for 

the  projectile  fired  at  = 4 seconds  is  about  20  ft  and  from  Fig.  4. 4-1  the 

go 

range  at  intercept  is  about  3600  ft,  the  range  to  the  target  when  the  projectile 
is  fired  is  5200  ft,  and  the  total  prediction  time  is  1.3  sec.  If  the  gun  pointing 
error  is  small  compared  to  the  predicted  position  error  and  if  the  target  does 
not  change  its  acceleration  drastically  during  the  prediction  time,  the  rms 
error  of  20  ft  found  here  is  a good  indication  of  the  miss  distance.  Knowing 
the  projectile  lethal  radius,  the  rms  miss  distance,  and  the  number  of  firings, 
the  kill  probability  can  be  estimated. 

Comparing  the  three  parts  of  Fig.  4.4-2,  it  is  clear  that  the  prediction 
error  is  very  sensitive  to  the  level  of  evasive  maneuver.  For  the  highly  evasive 
case,  the  firing  interval  during  which  the  error  is  less  than  100  ft  is  ai^ound 
four  seconds  for  Filters  B and  CE.  There  is  very  little  time  to  fire,  the  asso- 
ciated errors  are  relatively  large,  and  the  kill  probability  would  be  much 
lower  than  for  the  nonevasive  case. 

The  simulation  results  indicate  that  the  best  tracking  filter  for  these 
conditions  is  the  adaptive  b;uidwidtii  filter  (B).  It  is  clearly  best  for  the  non- 
evasive case  and  works  almost  as  well  as  any  of  the  others  at  higher  target 
maneuver  levels,  especially  when  rms  errors  under  100  fi  are  considered  to 
be  critical.  The  differences  between  ilie  performance  of  the  various  designs 
arc  not  great  compared  to  their  sensitivity  to  maneuver  level.  Depending  upon 
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the  lethal  radius  of  the  projectile  and  the  pohitinp:  accuracy  of,  Ihe  jjun,  these 
differences  may  or  may  not  have  a significant  effect  on  ovt^rall  kill  probability. 
In  addition,  large  prediction  errors  caused  by  the  target  changing  its  accelera- 
tion during  the  prediction  interval  (See  Fig.  4.3-2(c))  could  decrease  the  im- 
portance of  the  performance  differences  shown.  It  is  important  to  note  that 
much  better  accuracy  would  be  achieved  for  the  nonevasive  trajectory  if  the 
filters  were  designed  for  this  case  rather  than  for  the  nominal  trajectory;  this 
point  is  discussed  in  detail  in  Section  4.2.1. 


4.4.2  Sensitivity  to  Measurement  Noise  Level 

The  target  rms  predicted  position  errors  achieved  for  three  different 
noise  levels,  based  on  the  rms  estimation  errors  obtained  in  Section  4. 2.  z are 
given  in  Fig.  4.4-3.  The  rms  predicted  position  error  increases  significantly 
as  the  measurement  noise  increases,  as  would  be  expected.  The  choice  of  a 
’’best**  tracking  filter  for  this  set  of  conditions  is  not  clear.  Filters  A and  B 
are  superior  for  the  high  noise  case;  however  Filter  CE  is  best  in  the  low 
noise  case,  but  worst  in  the  high  noise  case.  These  results  taken  together 
with  those  presented  in  Fig.  4. 2-2  indicate  that  filter  B (adaptive  bandwidth) 
would  probably  be  the  best  choice  over  the  conditions  tested.  However,  in 
practice  the  best  choice  also  depends  upon  the  likelihood  of  the  off-nomlnal 
tracking  conditions  occurring. 


4.4.3  Sensitivity  to  Cor  re  1 alcd  Measurement  Noise 

Figure  4.4-4  shows  the  sensitivity  of  Filler  A to  correlated  measure- 
incnl  noise  for  the  numiiuil  tracking  situation,  using  the  results  obtained  in 
Section  4.2.3.  The  figure  shows  the  decrease  in  accuracy  that  (jccurs  if  the 
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(a)  rms  Measurement  Noise  Level  = 0.47  mrad  (Low  Noise  Case) 


(b)  rms  Measurement  Noise  Level  1.40  mrad  (Nominal  Noise  Case) 


Filmic  4.4-3  Prediction  Fn*or  Sensitivity  to  rms 
Measurement  Noise  l^*vel:  Nominal 
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EFijjure  4,4-4  Prediction  Error  Sensitivity  to  Measurement 

Noise  Correlation  for  Filter  A:  Nominal 
Trajectory 

D 

rms  measurement  noise  level  remains  constant  but  Us  correlation  time,  r,  is 
|j  changed  from  r = 0 to  r = 0.5  and  t = 1 sec^  Note  that  the  error  increases  by 

as  much  as  50%  over  the  range  of  r. 


4.4.4  Sensitivity  to  Delta  Rate 

Fiirurc  4.4-5  shows  the  improvement  in  prediction  accuracy  that  is 
obtained  tor  Filter  A in  the  nominal  tracking  situation  if  the  measurement  data 
rate  is  increased  fmm  10  measurements  per  second  to  20  and  40  measurements 
per  second.  Evidently  significant  improvement  is  possible  at  a higher  rate, 
but  whether  this  improvenienl  justifies  the  additional  computer  time  re- 
quired depends  u|)on  the  tradeoff  ijclween  increased  cost  and  improved  kill 
probability. 
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Figure  4, 4-5 
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Prediction  Error  Sensitivity  to  Measurement 
Data  Rate  for  Filter  A:  Nominal  Trajectory 


4.4.5  Sensllivliy  to  Target  and  Projectile  Velocity 


It  is  instniclive  to  see  how  the  prediction  results  change  if  different 
\.i*a.os  arc  chosen  for  the  projectile  velocity,  Vp,  and  the  target  closing  velocity, 
V.J  Figure  4.4-6  shows  the  results  for  Filter  A operating  in  the  nominal  track- 
inf*.  :’ifuaiion  where  v^  = 1000  ft/sec  and  Vp  is  varied  from  2000  to  3500  ft/sec. 
An  increased  projectile  velocity  helps  by  directly  reducing  the  required  predic- 
n tinu*  r,  *"  •*  rdven  time-to-go  and  hence  the  prediction  error.  Figure  4.4-7 
s»v»\7s  the  effect  of  ^^ariations  in  target  closing  velocity.  Increasing  v^  increases 
Ih*  prodiefii  Fj  ormr  for  a given  value  of  time-to-go  because  the  corresponding 
vFbie  of  tai/of  range  at  the  firing  instant  is  increased  (sec  Eq.  (4.3-5)).  There 
15*  a great  •ntrige  (o  targets  which  have  a very  high  closing  velocity,  as  well 
as  The  abilE'  • - maneuver  at  the  level  of  (he  nominal  Iraiicctorv  used  for  (his 
in%e.st  igatio:,. 
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Prediction  Error  Sensitivity  to  Projectile 
Velocity  for  Filter  A:  Nominal  Trajectory 


Figure  4.4-6 
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4.4.G  Summary  of  Prediction  Error  Results 

Several  general  conclusions  can  be  drawn  from  the  results  presented 
the  pretToding  sections.  First,  and  most  important,  tracking  filter  design 
I)'.  he  crucial  to  the  ultimate  effectiveness  of  the  gun  fire  control  system,  so 
cai'cful  and  thorough  design  effort  is  needed  for  each  specific  application. 

• 'sne  filler  design  of  either  the  fixed  or  adaptive  type  is  necessarily  best 
ill  target  engagement  situations.  The  fmal  selection  of  the  tracking 
't  M depends  very  much  on  the  complete  environmental  range  in  which  it 
operate,  the  speed  and  size  of  the  computer  available,  and  the  relative 
lance  assigned  to  tracking  accuracy  for  each  situation.  Even  for  a well 
• ified  set  of  trackuig  situations,  a single  filter  will  be  *^»e  best  for  each 
tt  L-r  of  the  set  and  the  final  design  must  be  select « on  the  ba^is  of 

t I V live  judgements  as  to  the  relative  importance  or  probability  of  each 
ion. 


Tracking  filler  comparisons  should  be  based  upon  prediction  error 
r than  estimation  error.  Certainly  these  performance  measures  are 
. 'o^ely  related,  but  prediction  is  ultimately  required  to  direct  the  gun  pointing 
and  prediction  error  depends  upon  prediction  time  and  the  target  motion 
•'  dfl,  as  well  as  the  tracking  filter’s  current  estimates  of  position,  velocity, 

•id  acceleration.  \Vliai  might  appear  to  be  a large  difference  in  performance 
i*‘  .u*en  two  fillers  based  upon  estimation  error  may  not  be  such  a large  difference 
ms  of  prediciion  error,  and  vice  versa,  depending  upon  the  attack  scenario. 

Tik  clioice  between  adaptive  and  fixed  designs  depends  ujx)n  the  spcci- 
M tuplic;?!  i.:fi.  Adaptive  dcsig7is  have  the  abili’y  to  offer  large  performance 
: ‘|)rov<  HI  • I:  m sfime  situations  but  thev  should  n(d  be  considered  as  the 
.\ei  b>  'I’l  pn>bloms.  Tliey  rciiiun*  additional  comj)le\il\  and  c(»m- 

p '»i  ion  ifMi-  and  iiia\  hi*  ovorlv  s»*nsilive  l*i  .aH'.\iaT»i*d  chnnia*s  in  the 
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innfturonii'>n(  noise  level.  Adaptive  traekini;  fillers  must  be  applied  cautiously 
and  stioulci  ahrays  be  judf;ed  against  a good  fixed  (nonadaptive)  design.  The 
raii'ic  of  possible  adaptive  methods  is  virtually  infinite  and  there  can  never  be 
cf-mjilete  assurance  that  the  best  desipi  for  a particular  application  has  been 
fou'id. 

The  results  in  this  chapter  indicate  that  guided  projectiles  are  needed 
t!.  gui'  fire  systems  which  are  required  to  have  a high  kill  probability  against 
t;.  with  a modest  manuver  capability.  Fipjre  4.4-20?)  shows  the  sive  of 
111*'  predicted  position  errors  for  the  various  tracking  filters  if  the  target's 
laas  H'  coloration  rate,  normal  to  the  line  of  sight,  is  approximately  1 g^sec. 

F o i.i.no-to-go  greater  than  six  seconds  all  of  the  filters  have  a rms  predic- 
. nor  greater  than  100  ft,  which  exceeds  the  lethal  radius  of  a conven- 
;i  r I projectile.  The  important  concept  here  is  the  following:  In  this  attack 
SC--C  ;rv),  the  target  will  be  verv  close  to  the  ship  before  the  rms  predicted 
i>o  or.in  error  is  smaller  than  the  lethal  radius  of  the  projectile.  This  means 
liwi  tlie  gun  is  effective  against  the  target  during  only  the  last  few  seconds  of  the 
ci.g..,rement.  Thus,  there  is  not  enough  time  to  fire  more  than  a few  limes  and 
evi  l,  il  the  target  is  hit  it  still  might  continue  on  its  trajectory  to  the  ship  and 
I serious  damage. 


There  is  an  ultimate  limit  on  the  prediction  ai  curacy  iwssilile  due  to 
the  llighi  lime  of  the  projcctilo.  Kvon  if  the  estinnles  from  the  tracking  filler 
ua-  perfc.l,  it  is  not  pos.sibtc  to  determine  the  targel'.s  exact  future  po.sition 
iiccause  lltc  target  can  maneuver  or  change  its  ;u  eclr  ration  during  tin  jirediction 


ioi'M  v,il.  As  iong  ,'ts  llic  target’s  future  lieliavior  is  un'iuiov.n  or  random  and  il 


a l.-Mgc  i-vasivc  capaliilil v.  il  i.s  impn<-si!>lv  lo  exactly  picdici  its  fuiuie 
|i Ill  this  case  the  mar.uilia!.  of  Ilie  miainien.  rms  pndicled  ixisiliim 
'•i:or  i„  iK-;ul\  imtopcmieiil  of  lln  l:-;-king  s-.stem  amt  (!i  pcn.ls  mosi  h Uiv.n  llu' 
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Beyond  these  general  conclusions,  the  results  show  that  the  specific 
adaptive  designs  tried  in  this  report  can  be  useful  for  gun  fire  control  systems. 
For  the  case  where  the  measurenient  noise  is  accurately  known,  the  adaptive 
bandwidth  and  hypothesis  testing  filters  give  significant  improvement  in  predic- 
tion error  over  the  fixed  design.  The  hypothesis  testing  design  is  especially 
useful  if  there  is  accurate  knowledge  of  the  relative  probability  of  different 
target  manuver  levels.  If  the  measurement  noise  is  poorly  knov..  and  can 
vary  as  much  as  a factor  of  three  above  and  below  the  nominal  design  level, 
the  advantage  of  the  adaptive  designs  is  partially  lost  because  they  are  more 
sensitive  to  measurement  nobe  level  than  the  fixed  designs.  However,  the 
adaptive  band  vidth  filter  still  performs  adequately  and  might  be  a good  design 
choice.  The  fixed  filter  design  also  has  some  good  features.  First  it  is  less 
complex  and  requires  much  less  computation  than  the  adaptive  designs. 

Second,  it  is  not  overly  sensitive  to  the  measurement  noise  level.  The  fixed 
filter  might  be  adequate  in  many  cases. 

Finally^  this  report  indicates  the  design  and  evaluation  procedure 
which  must  be  undertaken  to  select  a tracking  filter  for  a gun  fire  control 
system  application.  The  n siilis  show  quite  clearly  the  large  prediction  eri'ors 
that  can  result  if  the  filter  is  not  chosen  carefully.  Such  errors  can  seriously 
degrade  the  gun  system’s  kill  probability  and  must  be  kept  as  small  as  jxissible 
through  good  tracking  filter  design. 
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« Ftod  Filter  — Ttes  design  te  of  th»  target's 

Mtavlcr  aid  is  selected  to  be  tesnsttive  te  ct^^sss  ia 
that  beterior;  this  is  a gM^  desiei  *h»  the  amerces"  over 
the  range  of  tracking  sMsattem  aotic^adsd. 

• A<tolve  j^Bxta^h  Fiber  — This  davtga  adjoete  Ms  oeii 
bandwl^ith  bMed  ike  observed  target  betetvior;  t.e. , 
-■  it  tooetf  itself  on-line  to  aekieve  te^r^oved  trackteg 
accuracy  over  a range  of  Iraekfatg  sttmtkHis. 


ehleh  varions  fitter  dastgas  are  developed  aadfevahMtted. 

o . 

pEttiOBcldng  filter  for  this  problem  is  deserSMd  in  detaQ.  ^ rMSoe 
be  ia^teiMBted beeoaiM  clear;  l.e.,  in  c pndteal  teackhig  sibtt* 
generally  not  mough  known  aboiot  the  tax^  asd  tte  behavior  to 
. ^wW|tt  an  optoosB  tracldeg  Htter,  ftfeopthnam  decigns  with  bote  fliat::  a>*t 
additive  dttoetares  are  pn^osed  and  developed  as  practical  atterpatives  based 
vpoBtbekiaowtedgeavailriiileidbcmtMietergBC.  Ibear  dasigos  were  setected 
i^ter  a caretel  reviev  of  toe  technical  literatare  tor  tooeCmetooda  «eU  snked  _ 
to  preb||MD.  Tbi7  are: 
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Hypott^ste  Testing  Filter  — Tliis  design  ^ a bank  ol 
fillers  eadh  desl^d  for  a different  issumed 
level  <d  target  evasive  behavtor.  The  esttinatas  of 
each  filter  are  then  wei^ted  and  combined  to  give  an 
overall  estimate*  The  weights  are  ctetermined  from 
the  probabilities  that  each  of  the  filters  is  optimum* 

l^sid^  Testtog  Filter  — This  design  consists  of  the 
same  bank  of  fbmd  filters  used  in  the  hypothesis  test- 
ing filter  above*  The  final  estimate,  liewaver,  is 
selected  from  that  filter  \Kdiich  s^^ars  to  have  been 
operating  best  in  the  recent  past;  i.e. , that  filter 
whose  estimates  have  been  in  closest  agreement  with 
the  measured  data* 


A comparative  study  was  conducted  of  the  four  designs  described  above 
in  which  tracking  error  and  prediction  error  were  the  performance  measures*  - 
Tracking  error  was  defined  as  the  error  the  filter  makes  in  estimating  the  tir- 
get's  position,  velocity,  and  acceleration,  while  prediction  error  is  the  error 
tha  filter  makes  in  predicting  the  target’s  future  position*  The  advanta^a  and 
disadvantages  of  each  design  were  investigated  over  a wicfe^  range  of  i^tential 
tracking  situations*  The  sensitivity  of  estimation  and  prediction  accuracy  to 
changes  in  the  following  was  studied: 

• Target  evasive  maimver  level 

• Measurement  noise  level 

• Measurement  noise  correlation  time 

• Measurement  data  rate 

• Taiget  range 
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. is  8M«j^  to  be  cfmsL 
4,  i-4)  ia^dctes  IJKer' 

s^dfi^^l^ibetteb  tosn  fitter; 


AAptlv*  ftttsr  desigas  are  osoally  store 

the  smsareseeatBotoelerel  titan  find  deri^.  ti^ 

ctti  be  very  bxqtorttuit  if  the  actiea  level  li  cawsgl^ 


Fbr  almost  all  cases  Un^ed,  the  p:e^etlto  error 
rentte  todic^e  that  a gidded  pro^tlls  vmiM  be  ro* 
to  obtain  eeeepti^li  totimlity  aptJM  taititts 
that  manner  evastrely,  tatlen  toe  tarip^  to  to  be 
eUoeed  to  get  iritoto  a few  tboosmtd  f at  ^ toe 
before  tt  to  eedaged. 


Increased  measttremeat  rates  to^^it^  predh:^ 
accoracy  op  to  a point,  bat  restore  fawreased  rom*‘ 
ptdation  speed  in  tbe  tracktog  cmiqiator.  Ifbe  attl- 
mate  limit  is  prediction  aecnracy  Is  itoto?nrineif,tiot 

by  toe  tracktog  system,  bat  by  tbe  aMBtjr  o(  toe 

target  to  mamxver  darl  ig  the  fl^  tton  ef  Ob 
projecttfe. 


Ibe  8electo>n  trf  a tri^±‘3  Mter  reqato^* 
between  accarmnr  a^^^iaalty. 
dbUrns  reeatre  times  m^cam 

isKMae  oekipi  wMte  Itiey  aehdven  56^ 
re^ietloB  (at  maaener  levels)  in  toe  p(f« 
tion  error  for  soote  of  the  cases  investigated 
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TECHNOLOGICAL  FORECAST 


PROHIEM  DEFIMTION 


Tiu'  i)i  )bleni  addressed  in  this  report  is  that  of  designing  digital 
V king  ; l^oriUiins  for  shipboard  gun  fire  control  systems  which  can 
i raissile  targets  that  are  capable  of  unexpected,  large  evasive  maneu- 
. s.  Targets  having  the  capability  of  achieving  acceleration  and  accel- 

o *> 

. ♦ rate  levels  as  high  as  300  ft/sec^  and  150  ft/sec‘',  respectively, 

I . ' considered.  Th»  bjective  of  the  study  is  to  develop  tracking  algo- 
! )s  having  an  adaptive  structure,  permitting  automatic  tracking  of  a 
: . variety  of  target  maneuver  characteristics.  The  solutions  to  this 
. ! • .ioni  described  in  the  report  indicate  the  limitations  in  tracking  per- 
i .nuice  and  target  prediction  accuracy  achievable  with  modern  data  pro- 

i *’^g  technology.  This  information  will  be  valuable  to  the  Navy  for  de- 
» .jning  those  situations  in  which  hardware  improvements  — such  as 
accurate  tracking  sensors,  higher  measurement  data  rates,  and 
• ]>rojectiles  --  are  necessary  in  order  to  achieve  acceptable  projec- 
•.iss  distance. 


» > r.\TK  OF  TECHNOLOGY 

MiiC-  'rn  data  proces.sing  technology'  has  previously  been  applied  to 
i\  igiiiug  tnricct  tracking  algorithms  for  the  Navy’s  MK  80  Gun  Fire  Control 
.‘•x  ^lem.  Kalman  filtering  techniques  are  used  to  process  radar  tracking 
1 to  obtain  estimates  of  target  position,  velocity,  and  acceleration  for 
' in  pr(  dieting  the  target’s  future  position.  However,  the  target 


w 


maneuvers  investigated  in  those  applications  have  typically  been  at  levels 
on  the  order  of  one  ^ or  less.  The  emphasis  of  this  study  is  on  targets 
capable  of  much  higher  acceleration  levels,  and  capable  of  changing  their 
acceleration  quite  rapidly  in  an  evasive  fashion.  During  the  study,  adap- 
tive filtering  techniques  wore  developed  specifically  for  tracking  these 
types  of  targets  and  their  performance  was  evaluated  under  a wide  range 
of  tracking  conditions. 


6.3  SUGGESTIONS  AND  IMPLICATIONS 

Several  conclusions  obtained  from  this  study  arc  summarised  in  Sec- 
tion 5.2.  An  important  quantitative  I’esult  is  that  all  of  the  tracking  tccimiques 
yield  a predicted  root  mean  square  target  position  error  in  excess  of  50  feet 
for  evasive  targets  at  ranges  in  excess  of  5, 000  feet,  in  all  of  the  tracking 
situations  investigated.  This  is  near  the  lethal  radius  limit  for  conventional 
projectiles.  Some  improvement  in  prediction  accuracy  can  be  achieved  by 
employing  more  accurate  sensors  at  higher  ctita  rates;  the  former  may  re- 
quire the  use  of  infrared  or  optical  tracking  devices,  and  higher  data  rates 
will  require  increased  computer  caixibility  to  cycle  the  tracking  algorithm 
computalions  at  a faster  rate.  Alternatively,  higher  prediction  errors  can 
be  tolerated  if  guided  projectiles  with  the  caixibility  for  homing  on  the  tar- 
get are  developed.  The  output  of  this  study  includes  specification  of  algo-, 
rithm  design  techniques  needed  fur  tracking  maneuvering  targets,  and  in- 
dicates the  l irget  prediction  accuracy  limitations,  thus  providing  a basis 
for  judging  whether  hardware,  as  well  as  software,  imiirovements  are  re- 
quired to  acliicve  an  acceptable  probibility  of  kill  against  evading,  maneu- 
vering tar'':ets. 


As  stated  previously,  the  work  undertalien  here  wa.s  motivated  by 
the  increasing  threat  of  the  use  of  missiles  against  Na\7  shipping,  with 
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‘licir  aiu,iidant  hi^h  velociiics  and  aceelcrrilion  caiiabiliiies.  Our  effort  was 
t..  V rofi  toward  answering  the  question;  How  should  target  tracking  algo- 
•ihiiir  1.  designed  to  cope  with  maneuverable  targets  which  employ  unpi \ - 
Mijie  eviisive  actions  limited  only  by  their  maximum  acceleration  capa- 
i.  ..ty.  and  v;hat  arc  the  target  prediction  accuracy  limitations?  However, 

. ^ I no  ucar  term,  we  expect  that  the  guidance  laws  for  most  missiles  will 
i.coigtvjd  to  intercept  a ship,  without  having  a built-in  evasive  capability. 

‘ i sivucturc  of  these  guidance  laws  may  be  known  sufficiently  well  so  that 
' rr*  iiliig  performance  achieved  wHh  the  optimal  design  described  in 
or  3 can  be  closely  approximated  by  including  a mathematical  model 
•>.  ..iis5:ilc  guidance  law  v.’ithin  the  design  model  of  the  tracking  filter. 

• nvesiigatioii  of  the  target  prediction  accm'acy  achievable  under  these 
.amstanccs  merits  future  study. 
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APPENDIX  A 

OPTIMAL  FILTEmiNG  AND  PREDICTION  THEORY 


The  equations  of  the  opUmai  (KaP  .ivn)  illter  fcr  processing  noisy 
.'i:t  ahu*  tinent  data  are  reviewed  briefly  in  this  appendix  lo  provide  the  appro- 
pi  nle  t»:ickgroimd  and  notational  conventions  necessar>^  for  the  development  of 
.idaplive  tracking  filters  investigated  in  this  report.  Prediction  theory  is 
..5^  reviewed  and  its  relation  to  the  Kalman  filter  is  discussed.  The  devcl- 
ent  assumes  a basic  familiririly  with  random  variables  and  state  space 
.. . ,!inii,  for  linear  systems;  additional  detail  can  be  found  in  Ref.  18  and 
i 'Mpters  1 through  4 of  Ref.  19. 


A ! KALMAN  FILTER  EQUATIONS 

To  apply  Kalman  filtering  to  any  estimation  problem,  it  is  necessary 
•erive  a linear  stoc  hastic  first-order  vector  matrix  differential  equation 
.’  ici!  models  (he  maimer  in  which  the  system  states  interact  and  propagate  as 
,1  buK  lion  of  time.  This  equation  has  the  general  form 

x(t)  = F(t)x(t)  f G(t)u(t)  4 b(t)  (A.1-1) 

‘.•M*  \ ([)  is  an  m y 1 column  vector  representing  the  sy.stem  state,  F(l)  is  an 
m dynamics  matrix  which  defines  the  intoniction  of  the  state  vector  com- 
i.is,  and  u (1)  is  a p v 1 column  vector  of  white  gaussian  noise  inputs  such 

; « 

K|u(t)]  = 0:  Covfu)  ^ E|ii(I)u(t)^J  - Q(06(t-T) 

symbol  E [ ) doiu4os  niathemiilical  expectation:  Cov  (u)  denotes  the 
' •>vari:ince  matrix  of  ii . 
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The  matrix  G(t)  is  an  m v p distribution  matrix  which  indicates  how  each 
component  of  u(t)  affects  each  component  of  the  system  state  derivative,  x(t), 
and  b (t)  is  a r.i  y 1 column  vector  of  knowTi  system  inputs. 

In  the  target  tracking  problem  the  components  of  the  state  vector  x 
will  topically  include  position  anf  velocity  variables  plus  additional  components; 
representing  target  acceleration  and  perhaps  correlated  measurement  noise. 
Note  that  the  F,  G and  Q matrices  may  be  time- vailing.  In  the  target  tracking 
problem  the  system  state  equation  is  nonlinear,  unliiie  Eq.  (A.  1-1);  however, 
the  nonlincax*ities  can  usually  be  linearized  about  a nominal  trajectory  so  that 
the  linear  estimation  techniques  described  here  con  be  applied. 

The  solution  to  Fq.  (A.  1-1)  can  be  written  for  1 2;  in  the  form 
x(t)  = + J $(t,T)fG(T)u(r)  ^b(T)]d,•  (A.  1-2) 


where  x(tQ)  is  the  initial  value  of  the  system  state  vector  at  time  tg  and  ^(t,tQ). 
the  state  transition  matrix,  satisfies  the  matrix  differential  equation 


4.(i,  t(j)  = F(l)<h(t,  tg);  «(tj),  Ig)  = I 
When  F is  constant,  $ becomes  the  matrix  exponential. 


(A.  1-3) 


4>(l,  tp)  r e 


f«-V 


fA.I-4? 


In  a digital  tracking  system  one  is  generally  interested  in  the  stale  vector  at 
discrete  instants  of  lime.  Equation  (A.  1-2)  may  be  used  !o  relate  the  states  at 
two  instants  ^^f  lime,  t^^^  and  The  resulting  difference  equation  can  be 
written  in  the  form 


— n n-1  -n-1  — n-1  -n-i 


(A.  1-5) 
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where 


s = x((  ) 
-n  - n 


‘^n-1  - 'n-1^ 


f " 

^'n-1  " J t)G(t)u  (t)  d r 


n-1 


-n-1  = f 


(A.  1-6) 


n-1 


Observe  tlial  is  a naussian  wliitc  sequence.  That  is 


Efw^J  --  0 


E[w„ivy  = fO]:  j ^ 0 
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T 

1 w - 
-n-1  -n-1 
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Vr 


)Q(!)G(t)^  (t^.  r)dr 


n-1 


(A. 1-7) 


At  discrete  instants  of  lime,  measnremcnl.s  of  linear  combination.-;  of 
the  slate  variables  are  made.  The  equation  dcscribinii  this  measurortenl  j>rn- 
cess  has  (he  K^neral  form 


~ H X ♦ V 
-n  n -II  — n 


(A.l-«) 
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I S, / is  a vector  of  r raoasiired  quantiiio;;  at  lime  ( , II  i;:  an  r a m 
II  n n 

I V :r.i  r.iatrix  doscribiiu:  the  liiiv'ar  combinations  of  stale  variables  which 

. .•  ^ z in  the  absence  of  noise,  and  v is  an  r vector  of  zero  mean 
~n  -n 

;ara.sjicmcnl  errors  with  a covariance  matrix.  R at  time  t , defined 

an 
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At  any  time  t . !lie  objective  of  oplinml  estimation  liieo:^  is  to  process 
n 

' **,.  .jasiirtniiciiis  taken  up  to  that  time  to  produce  an  estimate  of  the 
.1.  ? tatc  havin*;  mininuuii  error,  in  a statistical  sense.  The  optimiza- 

<i  t • wrion  most  often  chg.scn  is  that  of  minimizing  the  mean  square  estima- 
t ..  • t . This  estimate  is  calculated  with  the  Kalman  filtcrlni^  aH;orithm. 

As  new  mc.asuremcnts  become  available  there  is  essentially  an  in- 

•u:;  change  in  onr  lmcwieflj!:e  of  the  state  x^.  Henoting  the  optimum 

• * * • of  X just  prior  to  the  a\*ai!abilltv  of  z as  x (-)  and  the  optimum 

— n * — n — n 

: ••  if  the  state  vector  immediately  after  processing:  is  b*e 

oo  Hltor  »;cneralcs  the  updated  optimum  estimate  of  the  system  slate 
< I *j  . ' > the  following  algorithm:^ 
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= ^-1  5n-l  ‘ :^n-r  J 


X fO  • .X  (-)  ♦ K - II  X (-1] 
-n  -n  n ' -n  ii  -n  ' 


(A.l-IOj 


Onlv  II-'  ili'^'-vctc  or  sninpicci  nicnsuroim-nt  form  ol  the  K.ilniun  filler  is  eoii- 
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where 


Xr  ^(c)  • the  optimal  estimate  of  ^ given  measuremei 


Zj.  Zg. 


X (-)  = the  optirn?!  estimrdc  of  x given  measurements 
-n  ^ -n 


Zj,  Zg.  ...  2„,J 


K = an  m^'  r ‘pin  matrix  furtlicr  descriV^ed  below 

n 


Observe  that  Eq.  {A.  1-9)  is  merely  the  discrete  form  of  the  contintoiis  tin;*' 
equal  ion 


x(t)  - F(t)  x(l)  I b(t) 


The  v]uantity  is  the  Kalman  gain  matrix.  Let  denote  liie  error 


made  in  estiniati  .g  x , i.c., 


X - X - X 

~n  -n  -n 


and  let 


P = Cov  X 
n -II 


^ rp  ”■ 

t-  I ''  ~ t 

- Ex  X 
[-n  -n 


Then  K is  computed  with  the  following  rccnrsi«>n  relutio.ns: 
n * 


K p (t)ir‘  ir’  = p (-)iE  m »>  (-ni*^ . n 

n n n n ii  r L n ii  •>  n J 


(A.l-M) 
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p (-)  - (j.  , p'  (t)  .:>  * f Q *, 

n ■ n-1  n-1  n-J  n-1 


p (4)  : p (-1  - P (-)•!'’  [}’,  1>  (-)H^  I R HP  (-) 

11  n n n n n n n n n 


a K H ) P {-)  (I  - K II  y 1 K H 

n 11  n n 11  n n n 


where  (he  (i)  and  (-)  noiiition  a*;:un  rel'eris  to  ll'e  quiaiiUies  before  and  alter 
measure^;:e’ii . 

Tne  ICn-’  ra'*  alf'orilhni  i.:a.  two  di.sU;ict  phase';.  Rqeutions  (A.  1-9)  and 
(A.  1-15)  du'ci  i’  e the  time  e.olution  of  fi.e  et^lc  e:  tlmato  and  itb  error  statis- 
tics lu'lwoen  a)cri.‘  uror.u  111:,  end  r l!ic  iaPueiice  of  system  dynamics  and  anise. 
This  piece.*  " is  eominonly  ieferred  * > uS  e.iix.TiOlr.tion.  Equations  (A.  1-tO) 
and  (A.  1-1C)  indicate  how  the  estimate  and  its  erro;  covariance  are 
at  the  measurement  time  to  reflect  the  now  iniorrnation  available.  The  extra- 
polation and  update  phases  of  the  Kalman  filter  are  summarized  in  Fi^.  A.  1-1. 
Information  flow  in  a lypieal  Kalman  filter  application  is  sliown  in  Fi^a  A.  1-2. 


Perhaps  the  most  uniqu  • feraure  of  the  Kalman  filter  is  that  ihe  per- 
formance anaiy.sis  of  the  filter  is  inherent  In  the  algorithm  for  K^.  The  niatri.x 
is  a eoiirplcle  de.se rip! ion  of  the  second-order  error  statistics.  Ii;  particu- 
lar, the  diagonal  terms  of  P^  rcprestMil  the  minimum  moan -.square  error  ob- 
tained in  cstimatin::  ear’i  comivueat  of  .\  . Note  thai  P is  specified  for  aU 

^ - n n . 

times  bv  E(is.  (A.  1-15)  and  (A.l-ldt.  Knowlodae  of  neiUier  x , x , nor  z 

-u  -n  ■ n 

is  required  to  obtain  a perfonnenee  analy.sis  for  (lie  up:  iia.il  filter. 

In  summary  the  following  c<nidi: ions  nui.''n  b*e  met  to  impleim  ai  an 


optiir.um  Kalman  filter: 


EXTRAPOLATION 


Pn(-)*‘I»n-1  *On-l 


UPDATE 

K„  - P„(-)hI;  [h„P„(-)hJ*R, 
p„H  = [i-k„h„]p„(-) 


THE  ANALYTIC  SCIENCES  CORP’ORATION 


• The  system  must  be  sufficiently  described  as  bein?: 
linear  and  obe\ir."  Eq.  (A.  1-1). 

• F(t),  G(t),  and  b (t)  must  be  knov/n  functions  of  time. 

• u (L)  must  be  a vector  of  zero  mean  t^aussian  white 
noise  inputs  with  known  covariance  ?natrix,  Q (t). 

e The  measurements  must  obey  Eq.  (A.  1-8)  and  H 
must  be  kno^\^l  for  all  n. 

• must  be  a vector  of  i^aussian  white  noise 
measurement  errors  witli  its  covariance  matrix, 

R^,  and  its  mean  knov/n. 

• To  initialize  the  filter  equations,  Eqs.  (A.  1-9)  and 
(A.  1-15),  Xq(-)  = E (Xq)  and  Pq(-)  = E fx^  x^  must 
be  provided. 

If  the  tracking  problem  met  all  of  the  conditions  above,  then  the 
design  of  an  optimum  tracking  filter  would  jusi  be  the  direct  implementation 
of  the  Kalman  filter  equations.  However,  in  practical  situations  one  or  more 
of  these  conditions  are  usually  violated  and  the  design  problem  is  more  dilfi- 
cult.  The  goal  is  then  to  desitpi  a suboptimal  filter,  based  upon  optimal  filter 
theory,  which  gives  nearly  optimal  performance.  Suboptimal  filtp»'  design  and 
evaluation  for  the  target  tracking  problem  are  the  major  topics  of  this  report. 


I 

A.  2 OPTIMAL  PREDICTION  THEORY 

‘ Given  eslimalcs  of  tlie  current  system  slate  it  is  ofte'i  required  Ihat 

I a prediction  bo  made  of  its  future  state.  In  Jnc  fire  control  problem  Ihe 

iUure  position  of  the  target  is  required  so  Uiat  the  projectile  will  intersect  the 


target  at  some  future  time.  To  aim  the  gun  properly  the  future  position  of  the 
target  must  be  known.  Kal.nan  filter  tlieorv  applies  directly  to  tliis  problem. 


THE  ANALYTIC  SCIENCE'S  CORPORATION 


{ ' 


I 


One  way  to  think  of  the  prediction  problem  is  to  look  at  how  the 
Kalman  filter  operates  wlion  no  measurements  are  being  made.  This  is  done 
by  setting  the  11^  matrix  in  the  measurement  equation,  Eq.  (A.  1-8),  equal  to 
zero;  then  Eq.  (A.  1-14)  indica^^s  "hat  the  Kalman  gain  matrix,  K^,  is  zero 
and  the  filter  equations  become 


X (-)  ^ X -(-)  4 b - 

-ir  ' n-.  n-r  -n-1 


(A.  2-1) 


P (-)  = , P ,(-)  4)^  , + 

n n-1  n-1  n-1 


(A.  2-2) 


Eq.  (A.  2-1)  is  the  optimum  prediction  algorithm  and  Eq.  (A.  2-2)  indicates  the 
propagation  of  the  prediction  error  covariance.  To  predict  the  future  system 


state  for  any  time  t^  beyond  the  present  time  these  equations  are  solved  itera- 


tively starting  with  the  current  state  estimate  and  its  covariance  matrix.  Note 
that  these  equations  are  independent  and  the  predicted  future  state  may  be  found 
without  using  Eq.  (A,  2-2), 


A-10 
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APPENDIX  B 

ADAPllVE  FILTERING  THEORY 

The  (lesii^a  of  an  optimal  Kalniiin  filter  to  estimate  the  states  of  a 
dynamic  system  requires  an  accurate  statistical  model  of  the  system  as  out- 
lined in  Appendix  A.  1.  This  model  completely  determines  the  filter’s  struc- 
ture vdiich  is  independent  of  the  measurement  data.  In  practice,  such  a model 
is  rarely  av^ailable  because  the  physical  system  is  not  well  knov/n  in  advance 
and  may  be  changirq^  unpredictably. 

One  approach  to  suboptimal  filter  design  is  to  specify  a fixed  filter 
design  and  e:q'»crimentally  determine  how  well  it  works  on  the  system  in  ques- 
tion. This  technique  is  satisfactory  in  many  situations  where  the  assumed 

system  model  is  a good  approximation  to  the  true  system  if  the  filter  is  do-  ^ 

■] 

signed  to  be  insensitive  to  errors  in  the  model.  However,  in  many  situations 
wheie  an  accurate  system  model  is  not  known  or  the  system  is  changing  in  an  i 

jj 

unkno^vn  manner  no  one  fixed  design  gives  adequate  performance.  It  is  this 
problem  which  motivates  the  search  for  adaptive  filter  structures  to  give  im-  j 

proved  performance. 

1. 

Since  tlie  measurement  data  contains  information  abort  the  system's 
structure  as  well  as  the  system's  state,  it  seems  reasonable  to  try  to  design  a 
filler  which  uses  the  data  to  identify  the  correct  system  model.  Essentially,  | 

an  adaptive  filter  is  a suboptimal  filter  with  a time-varying  structure  dependent  J 
upon  the  nicasurcmerl  data.  If  it  is  des*  .o<.*d  pro];crly,  it  sho\ild  be  able  to 
adjust  it:. elf  to  a previously  unknovm  or  chancing  sy.stcm  so  that  it  gives  near 
optimal  estimates.  The  degree  with  which  the  performance  of  an  adaptive 
filler  appi'(»achcs  that  of  an  optimal  filter  Ci.^pends  upon  how  fast  the  adaptation 

B-1 
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occurs  and  low  vrol!  Ihc  slnictnro  of  Hie  adaptive  filter  cun  model  the  structure 
of  the  optimal  filter.  Unfortunately,  there  is  often  no  unique  v/ay  to  select  a 
particular  adaptive  slructiu'c  for  a pariicular  application.  The  literature 
describes  many  different  types  of  adaptive  schorr.cs  and  the  only  way  to  judjz;e 
their  performance  is  by  experimentally  applying;  each  to  the  design  problem 
at  hand.  This  appendix  gives  the  details  of  two  of  the  three  adaptive  filtering 
techniques  tested  in  this  report.  The  third,  the  adaptive  bandwid"''  filter,  is 
discussed  in  Section  3.3.1. 

B.  1 HYPOTHESIS  TESTING  FILTER 

The  desi‘ai  philosophy  of  the  hypothesis  testing  filter  is  to  operate  a 
number  of  different  filters  in  parallel  and  compare  their  estimates.  This 
comparison  is  used  to  calculate  a final  estimate  which  is  a weighted  average 
of  the  individual  estimates;  those  which  appear  to  be  the  most  accurate  are 
given  the  largest  weights.  Specifically,  assume  that  over  some  time  interval, 


T,  the  optinnin  Kalman  filter  for  a particular  linear  system  is  one  of 
N kno^\^l  fillers.  If  this  assumption  is  correct,  then  it  is  possible  to  calculate 
the  probability  that  each  of  the  filters  is  optimum  over  that  interval.  The  state 
estimate  produced  by  the  hypothesis  testing  'Titer  is  the  sum  of  the  estimates 
from  each  of  Ihe  N individual  filters  weighted  by  the  probability  that  each  is 
optimum.  Tliis  desigm  is  theoretically  sound  and  can  be  expected  to  \wrk  very 
well  as  long  as  the  assumptions  upon  whicli  it  is  based  are  met.  The  primary 
disadvantage  of  (he  method  is  the  requirement  that  N Kalman  filters  be  imph^- 
mented  in  paralhd.  If  N is  large  and/or  the  order  of  each  filter  is  large,  this 
rcqiiiiemcnl  means  (hat  the  amount  of  computation  required  can  easily  become 
excessive.  For  the  purposes  of  this  report  a simplified  form  of  the  general 
hypothesis  testing  filter  is  developed  in  this  section.  A much  more  general 
discussion  can  be  found  in  Ref.  4. 
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Assume  iluit  the  system  to  v.hieli  tliis  method  is  apijlied  can  be  modeled 
by  Eq.  (A.  1-5)  with  the  deterministic  input  vector  b^  - 0,  that  is,  only  a vector 
of  "white”  -aussian  noise  sequences,  w^,  drives  the  system.  This  assumption 
represents  some  loss  of  !;eneralily,  but  many  systems  with  unknowi  inputs  can 
be  modeled  in  this  way.  If  b^^  is  knouai,  hov/ever,  it  can  easily  be  included 
tlie  equations  which  follow.  The  system  equation  then  becomes 


in 


^n  = Vl^n-l^'^n-1 


(B.1-1) 


where 


r T 

L^n-l'^n-l 

= Si-1 

(B.2-2) 

Further  assume  that  for  the  time  hiterval  under  consideration,  the  matrix 

is  constant  and  designated  by  Q^,  wliicli  represents  its  true  value  over  the 
interval,  T.  Assume  that  Q**  is  not  specifically  known,  but  it  is  known  to  be 
one  of  the  members  of  a set  of  N possibilities;  i,e.. 


Q*  <iq',  q'.  . . . q”j 


For  eacli  of  these  possibilities  there  is  a distinct  Kalman  filter  design 
defined  by  the  equations 


■ <'«-l  -n-l 


k / X 

X .(+) 


(B.1-3) 


-n 


-n 


k - 1,  2,  ...  , N (B.1^4) 


vk 


= Vi 


(B.l-5) 
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= P^(-)  fll  p’^(-)  + R 1 * fB.1-6) 

n n n L II  n n nj 

pV)  = [l  - H 1p’'(-)  (B.1-7) 

n'  ' L n nJ  n 

til 

where  the  superscript  k indicates  the  k filter  in  the  set  of  N filters,  each 

with  covariance  calculated  assuininj^  = Q^.  These  equations  also  assume 

that  f|),  H,  and  II  are  known  and  the  same  for  each  filter.  For  simplicity,  the 
1 1 ■ 

filter  initial  conditions,  Pq(-)  J^nd  Xq(-),  are  the  same  for  each  k.  This  com- 
pletely defines  tiie  structure  of  the  N Kalman  filters. 


The  adaptive  part  of  the  design  is  based  upon  calculating  the  proba- 
bilities of  N hypotheses  conditioned  on  the  measui  ements  taken  during  the 
interval.  The  first  !i>qDothesis,  designated  H , is  that  Q equals  Q . The 

2+2  t 

second  hypothesis,  H , is  that  Q = Q , and  so  on.  In  other  words,  li  is  the 

hypothesis  that  the  k^*^  filter  in  the  set  is  optimum  for  the  given  measurement 

data.  Once  the  probability  of  each  hypothesis  is  known,  the  optimum  state 

estimate  is  a linear  combmation  of  the  estimates  from  all  N individual  filters, 

calculated  by  multiplying  the  estimate  from  each  filter  by  the  probability  that 

the  hypothesis  corresponding  to  that  filter  is  true  and  then  adding  all  of  these 

weighted  estimates.  A diagram  of  this  structure  is  shown  in  Fig.  B.  1-1. 

To  discuss  how  the  various  probabilities  arc  calculated  it  is  necessary 
first  to  define  the  interval  measurement  history,  Z^.  Let  be  the  set  of  all 
measurements  taken  during  the  interval  in  question  up  to  and  including  the  m^*^ 
measurement,  taken  at  time  1^, 

n I -n-m-r  ’ -n-1  -n  t 

where  z _ - is  the  first  measurement  tahoii  in  the  interval  at  t 

n-m-j  ^ n-m-1 

The  probability  that  hypothesis  H is  true.  coiidilioii(?d  on  Ihc  measiircnicnl 
history  Z^,  is  designated  by  p | '/J  . 
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Q 0^'* 


Fip:urc  13.1-1  structure  of  Hypothesis  Testing  Filter 

Given  the  estimates  from  each  of  the  filters,  the  meiisurcmcnt  data, 
and  the  a priori  probability  of  each  hN^jothesis  at  the  beginning  of  the  intoi  /al, 

it  is  possible  recursively  calculate  p|ll^i  Z I . Using  Hayes’  rule  (Ref.  4.) 

th  * 

the  formula  for  (he  h proba!)ility  is 


(n.1-8) 


k 1 


where  f r-dimcnsional  normal  jirobability  density  func- 

tio’  r !)eiiig  llie  dimension  (»f  tlie  uu^asureinent  vector  The  value  of  tliis 
function  is  given  by 
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whc  re 


< = »n  ^ \ 


(B.1-10) 


and  [S^  I indicates  (he  determinant  of  the  matrix  sj|.  Note  that  H*'  is  the  k‘*' 
hypothesis  while  is  the  slate  obserTOtion  matrix.  In  this  manner,  p [h’'|  z *• 
IS  calcuLnted  for  each  k at  each  .samplLif^  time  in  (he  interval,  T.  "* 

The  overall  estimate  for  the  adaptive  filler  immediately  after  a 
nie.'isurcment  is  piven  by 


1:  = 1 


(B.1-11) 


with  its  corresponding  covariance  given  by 


N 

V*'-  Z)i-„‘wpKp„} 


k=  1 


(B.1-12) 


Up  to  this  point,  all  of  the  development  has  been  with  respect  (c  a 
single  interval  of  length  T.  Theort 'ically,  if  only  one  of  (ho  N fillers  is 
optimum  for  all  time,  T should  be  set  equal  to  infinity  and  the  design  is  com- 
plete. The  probability  of  the  correct  hypothesis  will  converge  to  one  :md  the 
final  estimate  will  be  equal  to  the  estimate  from  the  optimum  individual  filter 
If,  however,  the  correct  hypothesis  changes  with  time  because  the  correct 
model  for  the  system  changes,  then  a pmvision  inu.sl  be  made  for  the  adaptive 


D-6 


THC  ANALYTIC  SCfGNCLr»  COrtF'fDnATION 


u 


L. 


L 


r; 


r t 


I 


filter  to  restart  the  probability  calculation  to  test  for  a now  hypothesis  being 
true.  This  is  done  b>  iiiitiali/ing  the  fillers  rtncl  the  probiibility  calculation 
every  T seconds. 

At  the  beginning  of  the  first  interval,  the  N filters  must  be  initialized 
^ k k 

by  setting  each  Xq(-)  ^ priori  mean  and  covariance  of  the 

estimate.  That  is 


■>£< 


" [t  if] 


for  all  k 


for  all  k 


(B.1-13) 


(B.l-H) 


where 


k * k.  » 

X = X -X  (.) 
~0  ~0  ~0 


03.1-1 5) 


defines  the  initial  estimation  error  vector.  It  is  also  necossarv  to  choose  a set 
of  initial  values  for  the  pix>babilitics  calculated  recursively  using  Eq.  (B.l-R). 
This  choice  should  reflect  any  knowledge  available  about  the  system.  For 
example,  if  very  little  is  laiown  about  which  of  the  filters  might  be  optimum  it 
is  natural  to  choose  these  a priori  probabilities  equal.  After  the  first  interval 
and  at  the  end  of  each  interval  tlicrcafler,  the  filters  and  the  probability  calcu- 
lations are  reset  since  a new  h>qDothesis  might  hold  over  the  next  interval. 

Eacli  probabilily,  p |n’',  ] Is  read  to  Us  a priori  probability  airi  each  of 
the  N filters  is  re-initializcci  by  scUiit(;  its  estimate  equal  to  the  cojiihincd 
estimate  Riven  by  Eq.  (C.1-11)  and  by  scltin"  Us  eovarianee  nialrix  to  that 
Riven  l)y  Eq.  (H.  1-12).  The  entire  adaptive  filler  i.s  tlion  ready  to  c.alculatc  tlie 
probat)ilUy  of  each  Iiyjjothtsis  over  llie  next  interval  and  to  eontinuc  caleulatiiirt 
the  woiRhIed  estimate. 

The  dc.siRn  of  this  fil'or  requires  a number  of  chuice.s,  the  mo.st  im- 
portant of  wliich  i.s  the  set  of  hypotiiesos  or  Q*'  matrices.  TIh'  re.sel  interval, 
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the  a priori  probabililios  of  each  hypothesis,  and  (he  filler  iiiilial  condilions 
must  also  be  selecled  (o  oomplcfely  specify  the  dcsii^n.  These  choices  must 
reflect  what  is  knouni  in  adrrjice  about  (he  system  and  are  all  very  important  lo 
the  performance  of  the  ftlter. 

The  hypothesis  (eslinji  filter  has  some  primary  limitations  in  practice. 
First,  the  optimum  Kalman  filler  for  the  system  may  not  be  one  of  the  N possi- 
bilities over  a particular  interval.  This  could  be  due  lo  system  modeling; 
errors,  improper  choice  of  the  noise  covariance  matrix,  or  a choice  for 
the  set  that  does  not  cover  all  possibilities.  Secoiui,  the  choice  of  ihe 
reset  interval  ;uid  a priori  probabilities  is  often  i micwh.at  arbi  AU’y  because 
of  ignorance  of  the  syslem  dynamic  benavior.  Unforlimatcly  the  filter’s  per- 
formance is  oflcr  very  sensitive  to  lliese  effects.  In  particular,  as  discusse  ' 
in  Chapter  4,  if  ll.o  noise  covariance  matrix  is  incorrecl,  (he  adaptive 
al?,orithm  can  adapt  in  the  wron*;  directloti  and  f^ivr  (he  hipdie.si  weight  to  the 
estimates  from  that  individual  filter  which  is  producing  the  poorest  es(  .nates. 
The  final  difficully  with  this  desi^p)  is  the  requirement  t.hrd  N fillers  be  built  in 
parallel.  This  computational  burden  may  be  enacceplablo  in  many  applic.dior.s. 


B.  2 RESIDUAL  TESTING  FILTPJR 


The  inolivalioii  for  the  adaptive  design  described  in  this  section  is  Die 
sensitivity  of  the  hypothesis  testing  filler  lo  the  incorrect  choice  of  R^,  (ho 
mcasurenieid  noise  covariance  ir.alri:\  This  problem  is  n enlmncd  above  and 
demons! rated  draniatically  in  (he  rc.sulis  pre.sented  in  Section  4.2  and  4.4. 
Basically,  the  adaptive  algorithm  of  ihc  hypothesis  testing  filter  cepen  \ 
directly  upon,  R^^,  (see  Eep  (B.l-IO))  and  il  car  adapt  in  the  wnjiig  ”dircc(iun‘* 
if  R^^  is  m)(  t(pial  to  Mk  true  ludsc  covariance  matri.v.  'Fhe  res.dual  testing 
filter  has  Ihe  same  .sirarlure  w:i  hypothesis  les'uig  filler,  (see  Fig.  B.2-l)bul 
(he  adaptive  algorilltn;  i,-;  clior.un  cliffercnfly.  In  this  apprf),ach,  Uit  udat^livc 
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Fipirc  B. 2-1 


Structure  of  Bosic’ual  Tcsiin«i  Filter 


pari  is  (lcsi{’,ncd  so  that  it  is  irulcpcndent  of  the  clioicc  of  for  the  indivluual 
filters.  The  adaptive  al(;u.  llhiu  ohserves  all  of  (he  ebliniates  over  an  inierval 
T and  at  Die  end  of  the  interval  it  n.akes  a decision  as  to  which  filler  is  oporal* 
inicl>c*st.  Over  the  next  interval,  the  output  of  the  adaptive  al{;orit!nn  is  the 
output  of  ./.at  inciividiiul  filter  selected  as  best  over  the  previous  interval. 
Meanwhile  tlu?  (»I>servation  of  all  estimates  conlimies  so  that  a new  choice  can 
be  made  for  the  next  inierva' 


I ^ method  liy  which  the  choice  of  a best  fillei  is  made  at  tlr*  of 
cacti  intei  va  involves  Die  use  of  ific  residuaJ  sequonce  al  vixrh  filter,  the  ele- 
ments c)f  wliicli  are  defined  bv 


k " b I X 

V - II  X C-) 

•»*ii  -n  11  -n 


H-h 


(B.2-n 
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Tlii.*''  verior  ir>  prcniuK ipliefl  by  llic  KrJiv'rji  ^uin  lualrix,  k|^,  imcl  used  io  up- 
date ihc  state  e.stiiuate  aftci'  (he  moasuroment  is  made  (sec  Eq.  (A.  1-10)). 
Usin.a  t!ie  measuremeai  equeJ/ton,  Kq,  (A.  1-8),  and  Eq.  (A.  1-12),  the  equation 


fox'  V bccomoh 
~n 


~ H x'^(-)  + V 
-i\  n-i*.  “11 


(13.2-2) 


whore  estimation  error  of  the  filler  just  prior  to  r.rrival  of 

the  measuroment  z , A measure  of  each  filier’s  poiformancc  is  provided  bv 

k 

(ho  scalar  quantity,  defined  ns 


f'.i  = 21n 


T 

k*  . k 

Ay 
“11 


(B.2-3) 


whore  A is  my  syminetrie  positive  derinlte  maud::.  The  oypected  value  of  p 


denchod  by  is  [;iven  by 


Since  random  vnriaiiles  and  - 0,  it  follows 


(hat 


•-  !•:  II^  AH  x^(-)1  ! K [v  ‘'av  1 (B.2-4) 

‘n  L~'*  J L“’*  "'’J 


Krshiu'n  fillC’i'  111'"  'vv  .sUiit':!  tlial  ;;ii  (iiJhiia'i  fili'/r  miiiin.i^vs  !'  x *)>:; 

I '*  -fJ 

for  .nnv  jn»si(ive  scnii<Iefr;‘ile  mnlrix  H.  Not  ins  thai  the  matrix  1!  All  in  *•'(! 

11  II  * 

(n.2-*l)  is  |K>sitivc  stunide! inilr  because  A is  p.o  *ive  dellnili.  one  logical 
inenns  i»f  sclectiiq*  (!u  liest  lilJer  fL-  :i  Ml  of  N is  h*  select  (h:tl  wiiicii 
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minimizes  In  other  words,  since  an  opiimal  filter  is  known  fo  minin.izc 

rp 

E[x  ^ Bx]  and  the  first  term  on  the  right  side  of  Eq.  (B.2-4)  is  of  this  form, 

then  if  one  filter  out  of  a group  of  filters  is  found  to  have  a lower  than  all 

the  others,  it  is  a likely  choice  for  .he  optimum  filter.  If  none  of  the  filters 

— k 

is  optimum  then  that  filter  with  the  lowest  value  c>f  ^an  be  called  ’’best”. 

This  is  not  a rigorous  argument  for  choosing  tiie  best  filter,  but  it  is  a logical 

— k 

method  of  doing  so.  In  practice  is  not  !:nown  for  each  of  the  residual 

sequences,  but  an  estimate  of  it  can  be  made  over  some  interval  by  calculat- 

ing  the  time  average  of  g . Setting  A equal  to  the  identity  matrix,  the  esti- 
— k 

male  of  g is 
^11 


n +M-1 
o 


^k 


± 

M 


k 


T 


(B.2-5) 


n=n 


where  n^  is  the  value  of  lh(  time  index  n at  the  bcginnmg  of  an  interval  and  M 

is  the  total  number  of  sam])les  in  the  internal.  This  calculation  is  pcrfornied 

"~k 

for  each  of  the  N filters  and  the  filler  that  gives  the  minimum  g is  selected  as 
”besd*  over  the  next  interval. 

At  the  end  of  the  interval,  each  filter  is  rc-initializod  to  the  est  imate 
and  covariance  matrix  of  the  filter  which  was  judged  best.  Thus,  at  the  start 
of  each  i?itorval  all  of  the  N filters  have  the  same  initial  conditions.  This  is 
dor.e  to  reset  any  fillei  which  may  be  working  very  poorly  to  l!ie  estimate  and 
covariance  which  arc  probably  most  accurate. 


The  design  of  the  residual  testing  filter  consists  of  choosing  a set  of 
reasonable  filters  for  the  sy.stem,  initial  conditions  for  the  fillers,  and  a reset 
interval.  l iicso  clioices  must  be  made  based  upon  llie  particular  system  to 
which  the  liUor  is  applied  and  .should  reflect  any  kiiow’lcdge  available  abom  (lie 
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system . The  filter  has  essentially  tho  same  computational  requirements  as 
the  desipjfi  developed  in  the  last  section.  It  is  less  sensitive  to  an  incorrect 
choice  for  the  mep,suvement  noise  covariance;  i.e. , it  does  not  adapt  in  the 
wrong  direction.  However,  it  tends  to  be  less  accurate  than  the  hypothesis 
testing  filter  in  cases  where  the  design  pat'anielers  for  the  latter  are  accur- 
ately known. 
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